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Notation 

Cp  -  Pressure  coefficient  defined  in  eq.    (3) 

a  -  Yaw  angle  relative  to  laboratory  space 

<J)  -  Pitch  angle  relative  to  laboratory  space 

P.  -  Total   pressure 

P  -  Static  Pressure 

a  -  Yaw  angle  relative  to  probe  axis 


rp 


-  Pitch  angle  relative  to  probe  axis 


P       -  Pressure  indicated  by  the  probe 
a       -  Yaw  setting  of  the  probe 
<p        -  Pitch  setting  of  the  probe 


relating  to  probe  positions  shown 
II  in  fi qui 


Subscripts 

I 

II 

III  in  figure   3 
l 

IV  I 

A  relating  to  A  type  probe 

B  relating  to  B  type  probe 

Lo  lowest  value  in  an  array 

rp  relative  to   the  probe  axis 

Up  highest  value  in  an  array 


Notation  of  PENPTS 

DD         -  denominator  in   the  subsequent  line  of  the  program 

HT         -   penetration  height  (Cp  of  a  calibration  surface) 

IR         -   a  flag  with  a   value  of  either  1   or  2   to   identify  first 
or  second  penetration   point.     This   flag  also  controls 
return  of  values   to   calling  program  or  subroutine 

IB         -  a   flag  indicating  the  location  of  the  last  rectangle 
in  the  band  checked  for  a  penetration  point.     It  pre- 
vents  repetition  of  the  scanning  in  the  103  Do  loop 

NX         -   number  of  X  values   in  the  calibration  surface  grid  (ex  values  of  the 
calibration  surface) 

NY         -  number  of  Y  values  in  the   calibration  surface  grid  (     values  of  the 
calibration  surface) 

X(I)     -  X  values  in  the  grid 

XM        -  value  of  X  at  point  B   (fig.   5) 

XR(  I)  -   two  values  of  X  returned  by  the  subroutine  as  results 

XS         -  currently  calculated  value  of  X  which   is   the  present 
result  and  is  later  stored  in  XR( I ) 

Y(J)     -  Y  values   in  the  calibration  grid 

YC         -  intermediate  value  used  in  the  subsequent  line  of  the  program 

YG         -  the  Y  position  of  the  penetrating  line* 

Z(I.J)-  Values  of  the  Z  =  Cp  calibration  surface  above  grid  points 

ZX         -  value  of  Z  either  at  point  A  or  point  C   (see  fins.    5, 
6   and  7) 

ZM         -  value  of  Z  at  point  B 


*  It  should  be  pointea  out  here  that  the  penetrating  line 
is  parallel  to  the  X  axis  at  a  height  HT  =  Cp  above  the 
X,Y  plane  and  at  a  distance  YG  from  the  X  axis  on  the  X,Y  plane 


Notation  of  INTSCS 


AA       -  constant  defined  in  eq.  8 

AB       -  constant  defined  in  eq.  9 

ARIN      -  ordinate  at  which  scanning  starts 

ARLO      -  lowest  ordinate  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  all  calibration 
surfaces  used) 

ARM!      -  ordinate  measured  relative  to  probe  1  (the  probe 
belonging  to  the  first  curve  Cp  =  const.) 

ARM2      -  ordinate  measured  relative  to  probe  2  (the  probe 
belonging  to  the  second  curve  Cp  =  const) 

ARP1      -  setting  ordinate  of  probe  1 

ARP2      -  setting  ordinate  of  probe  2 

ARRES  (1)  -  abscissa  value  returned  to  INTSCS  from  PENPTS 

ARUP      -  highest  ordinate  of  calibration  surface  (as  a 

condition  it  must  be  identical  for  all  calibration 
surfaces  used) 

AR2R      -  ordinate  relative  to  laboratory  space 

AR2RJ     -  ordinate  relative  to  laboratory  space  of  previous 
scan 

BA       -  constant  defined  in  eq.  8 

BB       -  constant  defined  in  eq.  9 

DAR       -  ordinate  step  for  scanning 

HT1       -  penetration  height  of  calibration  surface  of  probe 

1  (Cp  of  probe  1 ) 

HT2       -  penetration  height  of  calibration  surface  of  probe 

2  (Cp  of  probe  2) 

IEPS  -     a  flag  which  is  equal    to  2  when  the  calculation  is 

carried  out  as  shown  in  fig.    11,  otherwise  it  is 
equal    to  1.      (This   information   is  required  in  the 
main   program  and  is   returned  to  it) 


ISL  -  a  flag  which   is  equal    to  10  if  one  of  the  left 

hand  side   results  returned  from  PENPTS  is  1000.0 
(i.e.,   is   not  a  penetration   point).     Otherwise 
it  is  equal    to  1 

ISR  -  a  flag  which  is   equal    to  10  if  one  of  the  riqht 

hand  side  results   returned  from  PENPTS  is   1000.0 
(i.e.,   is   not  a  penetration  point).     Otherwise 
it  is  equal    to  1 

K  -  a  flag  indicating  whether  the  first  or  the  second 

intersection  is  evaluated 

Nl  -  number  of  abscissas   in  calibration  surface  grid 

(fig.  4) 

N2  -  number  of  Ordinates  in  calibration  surface  grid 

(fig.  4) 

RES1(K)  -  evaluated  and  returned  abscissa  of  intersection 

RES2(K)  -  evaluated  and  returned  ordinate  of  intersection 

Rl  -  radius  in  fig.   11    for  probe  1 

R2  -  radius   in  fig.   11    for  probe  2 

RAIL  -  abscissa  of  point  AIL  in  fig.   9  or  10 

RAIR  -  abscissa  of  point  AIR  in  fig.  9  or  10 

RAJL  -  abscissa  of  point  AJL   in  fig.   9  or  10 

RAJR  -  abscissa  of  point  AJR  in  fig.  9  or  10 

RBIL  -  abscissa  of  point  BIL  in  fig.  9  or  10 

RBIR  -  abscissa  of  point  BIR  in  fig.  9  or  10 

RBJL  -  abscissa  of  point  BJL  in   fig.   9  or  10 

RBJR  -  abscissa  of  point  BJR  in  fig.   9  or  10 

XI (I)  -  abscissas  of  calibration  surface  of  probe  1   or 

its   transformation  as  determined  by  calling 
statement  in  main  program 

X2( I)  -  abscissas  of  calibration  surface  of  probe  2  or 

its   transformation  as   determined  by  calling 
statement  in  main  program 


X3(I) 

X4(I) 

Y1(I) 

Y2(I) 

Y3(I) 

Y4(I) 

Z1(I,J) 

Z2(I,J) 

Z3(I,J) 

Z4(I,J) 


-  abscissas  of  calibration  surface  of  probe  1   or  its 
transformation  as   determined  by  calling  statement 
in  main  program 

-  abscissas  of  calibration  surface  of  probe  2  or  its 
transformation  as   determined  by  calling  statement 
in  main  program 

-  ordi nates  matching  XI (I 

-  ordi nates  matching  X2(I 

-  ordi  nates  matching  X3(I 

-  ordi nates  matching  X4(I 


-  values  of  Cp  for  XI  (I 

-  values  of  Cp  for  X2( I 

-  values  of  Cp  for  X3(I 

-  values  of  Cp  for  X4( I 


and  Yl ( I) 

and  Y2(I) 
and  Y3(I) 
and  Y4(I) 


Notation  of  Main   Program  VDR 

A  -   variables   in  library  routine   IXCLOK  used   to   evaluate 

computation   times 

ALF         -   values  of     a     returned  to  VDR   from   INTSCS 

ALFA       -  a     input  when  experiment  simulation   is   carried  out 


ALFC 


-    a 


rpIII 


ALFCN     -   new  value  of    ^rDTjT   for  next  iteration 

ALFD      -«rPIV 

ALFDN     -   new  value  of     a     IV  for  next  iteration 


ALF1 
ALF2 
ALF3 
ALF4 


"  aI 


aII 


aI 


"  "III 

AM    -  mach  number 

CPA    -  Cprp, 
CPB    -  CprpH 

cpc   -  cPrpIII 

CPD    -  CprpIV 

EPSPS     -   relative  difference  between  present  static  pressure 
and  its  value  in  the  previous   iteration 

EPSPSG  -  convergence  criterion   (on   the  static  pressure) 

FALF       -  final     a     value   (at  convergence) 

FPHI       -  final      4    value  (at  convergence) 

IA50       -  a   flag  in   IXCLOK 

IA60       -  a  flag  in   IXCLOK 

ICOPS     -  number  of  PS   corrections  carried  out  to  ensure  the  Cp's 
are  in  calibration  range 

IEPS       -  a  flag  governing  the  convergence  criterion  returned 
from  INTSCS 


I  EPS!     -  a  flag  governing  the  convergence  criterion  for  probes 
I  and  II 

IEPS2     -  a  flag  governing  the  convergence  criterion  for  probes 
I   and  III 

I  IT        -  number  of  iterations  on   p 

ISCAN     -  number  of  static  pressure  scans  from  initial   static 
pressure  guess  upwards 

IXCLOK  -  system  subroutine  for  computation  time  evaluation 

NOCOPS  -  a  flag  which  when  equal    to  1   causes  static  pressure 
corrections  to  be  carried  out  such   that  Cp's  are 
always   inside  calibration   range,   and  when  equal    to  2 
cause   these  corrections   to  be  skipped 

NOSIM  -  a  flag  which  when  equal  to  1  causes  the  program  to 
simulate  velocity  measurement  experiments  and  when 
equal    to  2  causes   the   program  to   reduce  measured  data 

NX  -  number  of  calibration  surface  matrix  abscissa's, 

a  values   (must  be  identical    for  all    calibration 
surface  matrices  used) 

NY  -  number  of  calibration  surface  matrix  ordi nates, 

<t>  values   (must  be  identical    for  all    calibration  surface 
matrices  used) 

PDYN       -  dynamic  pressure 

PHI         -  value  of     <$>     returned  to  VDR  from  INTSCS 

PHIB       -    ^     TT 
TpII 

PHII       -    <f>  input  when  experiment  simulation  is   carried  out 
PHU        -    Vpl 


PHI2 
PHI3 


VpII 
Vpl 


PHI4     -  Vpin 

PPA         -   signal   of  probe  in  position   I 
PPB         -  signal   of  probe  in  position    II 


PPC  -  signal   of  probe  in  position  III 

PPC  -  signal   of  probe  in  position   IV 

PS  -  static  pressure 

PSIN  -  initial   guess  of  static  pressure 

PSN  -  new  value  of  static  pressure  for  next  iteration 

PSC  -  corrected  static  pressure 

PST  -  memorized  corrected  static  pressure 

PT  -  total    pressure  measured 

RELXPS  -  relaxation  factor  for  convergence  on  the  static 
pressure 

VI RTIM  -  virtual   computation  time 

XA  -  abscissas  of  calibration  surface  matrix  of  probe  A 

XAX        -  transformed  abscissas  of  calibration  surface  matrix 
of  probe  A 

XB  -  abscissas  of  calibration  surface  matrix  of  probe  B 

XBX        -  transformed  abscissas  of  calibration  surface  matrix 
of  probe  B 

XPC        -  yaw  setting  of  position  III 

XPD        -  yaw  setting  of  position   IV 

XRIN      -  a  value  at  which  scanning  starts 

XRLU       -  lowest     a     value  of  calibration  surface  (as  a  condition 
it  must  be  identical    for  probes  A  and  B) 

XRUP       -  highest     a     value  of  calibration  surface   (as  a 

condition  it  must  be  identical    for  probes  A  and  B) 

YA  -  ordi nates  of  calibration  surface  matrix  of  probe  A 

YAX         -   transformed  ordi nates  of  calibration  surface  matrix 
of  probe  A 

YB  -  ordi nates  of  calibration  surface  matrix  of  probe  B 

YBX         -  transformed  ordi  nates  of  calibration  surface  matrix 
of  probe  B 


YPA  -  pitch  setting  of  probe  A 

YPB  -  pitch  setting  of  probe  B 

YRIN  -   <p     value  at  which  scanning  starts 

YRLO  -  lowest     <j>    value  of  calibration  surface   (as  a 

condition   it  must  be  identical    for  probes  A  and  B) 

YRUP  -  highest     <f     value  of  calibration  surface   (as  a 

condition   it  must  be  identical    for  probes  A  and  B) 

ZA  -  Cp     .   values  of  probe  A 

ZAX  -   transformed  Cp     .   values  of  probe  A 

ZAMAX  -  maximum  of  ZA  array 

ZAMIN  -  minimum  of  ZA  array 

ZB  -  Cp     «  values  of  probe  B 

ZBX  -  transformed  Cp     R  values  of  probe  B 

ZBMAX  -  maximum  of  ZB  array 

ZBMIN  -  minimum  of  ZB  array 


Introduction 

Experimental  knowledge  of  the  flow  field  generated 
by  rotating  turbo  impellers  is  of  prime  importance  in  the 
research  and  development  of  turbomachi nery .   It  is  essential 
for  the  refinement  of  design  methods,  for  the  development 
of  new  flow  models  which  include  secondary  flow  and  tip 
clearance  effects,  and  particularly  for  the  verification 
of  new  computer  codes  developed  to  calculate  the  flow  through 
rotating  blade  rows. 

In  recent  years  laser  veloci meter  techniques  have 
been  applied  successfully  to  measure  the  flow  both  inside 
and  downstream  of  rotors.  (Ref.  1  for  example).,   It  has 
become  clear  however,  that  the  laser  techniques  are  only 
reliable  in  the  hands  of  experienced  investigators.   A 
window  which  remains  clean  is  essential,  and  seeding  is 
usually  required.   Laser  techniques  do  not  measure  the 
pressure  field  and  usually  can  only  measure  two  components 
of  the  velocity  unless  the  axis  of  the  laser  is  tilted. 
Difficulty  is  also  encountered  when  measuring  close  to 
walls.   Hence  there  are  reasons  to  consider  alternative 
techniques,  particularly  if  they  are  simpler  to  apply  routinely 
in  stationary  turbomachi nery  passages.   Furthermore,  the 
achievement  of  redundancy  in  measuring  the  flow  field 
behind  the  impeller  is  itself  a  worthwhile  goal.   The  present 
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work  deals  with  the  application  of  a  particular  system 
of  small  high  response  pressures  probes  at  the  exit  of 
an  impeller. 

Measurements  behind  an  impeller,  in  the  stationary 
bladeless  gap,  are  simpler  to  make  than  measurements  within 
the  rotating  passages.   Transducer  probes  can  be  installed 
through  the  stationary  machine  casing  and  the  data  transmitted 
without  resort  to  slip  rings  or  rotary  transformers.   The 
sensor  is  not  subjected  to  the  centrifugal  field  or  to 
the  vibration  of  the  rotor.   However,  the  flow  to  be  measured 
is  then  fluctuating  at  bl ade- passi ng  frequency  and  any 
system  of  sensors  must  be  calibrated  for  a  wide  range  of 
possible  Mach  number,  pitch  angle,  yaw  angle  and  pressure 
variation  -  and  yet  must  be  capable  of  the  necessary  fre- 
quency response. 

In  Ref.  2,  a  method  was  described  for  using  two  semi- 
conductor pressure  probes  together  with  the  technique  of 
synchronized  sampling,  to  obtain  the  distribution  of  the 
velocity  vector  downstream  of  a  rotor.   The  geometries 
of  the  two  probes,  designated  Type  A  and  Type  B,  and  their 
installation  in  the  compressor  annulus  are  shown  in  Fig 
1.   It  was  argued  that,  in  principle,  by  rotating  the  probes 
in  yaw  about  their  tips  and  controlling  the  sampling  of 
the  data  from  each  probe  to  be  at  the  same  position  in 
the  rotor  frame,  the  system  of  two  separate  probes  could 
be  used  to  acquire  data  at  a  point  in  the  periodic  flow 
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FIGURE    I.      TYPE   A   AND   TYPE    B    PROBES   AND   COMPRESSOR    INSTALLATION. 
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from    the    rotor,    corresponding    to    data    normally    obtained 
from    the    multiple    sensors    of    four-    or    five    hole    pneumatic 
probes    when    measuring    velocity    in    steady    flows.      The    tech- 
nique   promised    the    use    of    probes    having    the    simplest    geo- 
metry   and    thus    avoided    the    large    size,    expense    and    unreli- 
ability   of   multiple    sensor    probes    which    incorporate    multiple 
semi-conductor    transducers    (Ref.    3).      Because   of    the    simple 
sensor    tip    geometry    (that    of    a    cylinder    at    incidence    to 
the    flow)    the    unsteady    response    was    likely    to    be    as    good 
as    could    be    expected    of    any    single    physical    sensor. 

The    two-probe    technique    is    strictly    applicable    only 
to    periodic    flows.      However,    data    obtained   on    successive 
rotations    of    the    rotor    can    be    averaged    to    eliminate    fluctua- 
tions   which    are    not    periodic.      This   was    shown    to    be    effect- 
ive   in    tests    reported    in    Ref.    2    in    which    a    single    Type 
A    probe   was    used    to    establish    the    peripheral    blade-to-blade 
distribution    of    flow   yaw    angle. 

In    order    to    obtain    velocity    from    the    pressure    measure- 
ments   which    can    be    obtained    from    the    two    probes,    the    steady 
response    characteristics    must    first    be    established    in    cali- 
bration   tests    carried    out    in    a    known,    controlled,    uniform 
flow.      Second,    a    method    of    applying    the    calibration    to 
measurements    made    in    an    unknown    flow    must    be    devised. 
In    the    present    method    two    different    approaches    have    been 
followed.       In    Ref.    4,    a    technique    is    given    for    representing 
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and   applying    the    probe    calibration    analytically.    When    first 
applied    the    method    gave       surprisingly   good   accuracy    (1    -   2%)    since 
the  method    required   that    the   probes   had   characteristics 
which    could    be    well    represented    analytically.      This    in 
turn    required    that    the    probe    tips    be    geometrically    precise, 
a    feature   which   was    not    present    in    the    first    generation    of 
probes  . 

A    second,    more    general    approach    is    reported    here,    wherein 
the    calibration    of    each    probe    is    represented    by    a    two-dimen- 
sional   array   of   pressure    coefficients.      The    application    of 
the    calibration    given    in    this    form,    in    an    unknown    flow    required 
the    development   of    special    numerical    procedures.      The    purpose 
of   the    present    report    is    to    document    the    analysis    and    the 
Fortran    program    developed    to    apply    the    method. 

In    its    present    form,    the    method    does    not    require    that 
the    calibration    "Surfaces"    be    symmetrical    about    any    axis    or 
be    expressed    in    analytical    form,    but    does    require    that    the 
pressure    coefficient   be    independent   of    Mach    number.      The    lat- 
ter   restriction    could    undoubtedly    be    removed    by    introducing 
additional    iterative    steps.       Further,    in    the    present    method 
only    five    measurements    have    to    be    taken    to    determine    uniquely 
a    velocity    vector    at    a    point.      Throughout    the    report,    the 
Fortran    program    notation    has    been    used    to    describe    the    physics 
and   equations    involved    in    the    solution. 
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Mathematical    model 

Assume    the   A    and   B    type    probes    of    ref.    [1]    (see    also    fig.    1) 
to    be    immersed    in    a    three-dimensional    steady    flow    field.*      The 
pressure    response    of    each    of    these    probes    in    given    gas    is    func- 
tionally   described    by    four    variables    as: 


Pp   =    Pp    (a,<|>   ,    PT,    Ps) 


(1) 


If  a  pressure  coefficient  is  defined  as 


P  -  P 

_  _j> 1 

Cp     "  P,  -  P 

Hrp     T    s 


(2) 


The    calibration    surface    of    each    probe    is    given    in    the    general 
case    in    form   of    a    matrix    of    values    of    Cp,    where 


Cp         =    Cp       (a      ,  d>      ) 
prp  Krp\     rp  >  ypp/ 


(3) 


The  pressure  coefficient  defined  in  this  way  has  only  a  second 
order  dependence  on  Mach  and  Reynolds  numbers  in  the  range  of 


*  For  our  purpose,  using  the  synchronized  sampling,  the  flow 
field  behind  the  impeller  is  steady,  although  the  probes  re- 
quire a  high  speed  response  because  of  fluctuations. 
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0    <Ma<    0.7      in    turbulent    flows    (Ref.    2),    so    that,    to    first 
order    their    influence    is    neglected    in    writing    eq .    (3). 

If    the    type    A    probe    is    rotated    about    its    axis    into    three 
different    positions    (i    =    I,    III,    IV)    readings    are    taken,    and 
the    type    B    probe    is    fixed    in    position    II    and    a    single    reading 
is    taken,    the    following    four    equations    can    be    written. 


CP 


Ai 


pAi 


PT    -    pr  ^Ai 


=    C  p  fl  •  ( a  ,  cj> ) 


Cp 


P  P 

pBII         rS 


(4) 


BII  PT    -    P5 


CpB  j  j  (a,(f> ) 


it    should    be    pointed    out    here,    to    avoid    misunderstanding    that    a 
and     <J>      are    defined    in    a    coordinate    system    relative    to    the 
machine    axis    and    not    relative    to    the    probe    axis.      In    the    set 
of    four    equations    (4)    there    are    four    unknown    quantities,    namely: 
a,    <(>,    P-j-    and    P^    .      These    are    the    quantities    to    be    evaluated    using 
the    measured    data.      Together   with    the    stagnation    temperature    they 
define    the    flow    field    uniquely.      The    four    equations,    resulting 
from    the    four   measurements,    should    be    sufficient    to    determine 
the    four    unknown    quantities. 

But    the    problem    is    complicated    by    three    facts: 

1)  The    calibration    surfaces    are    not    generally    known    in 
analytical    form. 

2)  The    calibration    surfaces    are    double    valued    in      a      and      <f 
i.e.,    for   a    given    Cp    and      a      there    exist    two       <j> 
values,    or    for    a    given    Cp    and      41      there    are    two      a 
values    which    satisfy    eq.    3. 
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3)   As  a  result  of  the  measurement  the  value  of  P   rather  than 

of  C   is  determined. 
P 

Since  the  calibration  surfaces  are  not  given  in  a  simple  analytical 

form  the  solution  has  to  be  numerical.   An  iterative  procedure  is 

required  because  P   and  not  C   values  are  measured.   First  P   and 

P  P  s 

P-j.  have  to  be  guessed  to  yield  C   values,  knowing  the  four  measured 
values  of  P  ,  and  using  eq .  A.   The  guess  is  then  iteratively 
corrected  to  converge  on  the  solution.   However,  convergence  of  the 
iterative  procedure  is  complicated  because  of  the  double  valued 
nature  of  the  calibration  surfaces. 

This  method  of  iteration  shown  in  Fig.  2  was  attempted  initially 
for  the  evaluation  of  a  measured  point.   In  practice  the  initially 
suggested  procedure  converged  in  some  cases  and  diverged  in  others, 
depending  on  the  values,  and  signs  of  a  and  <J> .   This  was  not 
surprising  as  convergence  on  two  variables  is  not  likely  to  be  a 
simple  matter.   However,  in  order  that  the  measurement  technique 
be  useful,  convergence  on  the  correct  solution  for  a  general  set  of 
measurements,  is  absolutely  necessary.   In  practice,  this  can 
certainly  be  achieved  if  one  of  the  two  iteration  variables  is 
obtained  by  measurement.   Since  the  static  pressure  is  a  difficult 
quantity  to  measure  even  in  a  steady  flow  field,  only  the  stagna- 
tion pressure  measurement  need  be  considered.   It  is  possible  that 
the  time-varying  stagnation  pressure  could  be  measured  with  a  suit- 
ably designed  Kiel  probe.   Data  would  then  be  taken  by  synchronized 
sampling  from  the  fixed  Kiel  probe,  from  the  type  A  probe 
rotated  into  two  positions,  and  from  the 
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GUESS 

P     £  P 

T    *■     S 


COMPUTE 
FOUR   C     VALUES 
FROM    FOUR    Pp  VALUES 


EVALUATE 


CALCULATE 

P     $  P 
KT  q     S 


CHECK 
FOR    CONVERGENCE 
OF    PT  g    P$ 
SIMULTANEOUSLY 


Fig.    2:      Unstable    iterative    solution    using    measurement 
of   A    type    and    B    type    probes    only. 


fixed  type  B  probe.   The  method  of  solution  is  then  as  shown 

in  Fig.  3.   The  method  shown  in  Fig.  3  proved  to  converge  under  all 

conditions.   It  is  described  in  detail  in  the  following  pages. 
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GUESS 


COMPUTE 
4    Cp    VALUES    FROM   4 

Pp  VALUES    t  MEASURED  P 
VALUE 


EVALUATE 

a  g   $ 


CALCULATE 


CHECK 
FOR  CONVERGENCE 
ON    R. 


Fig 


Stable  iterative  solution  using  measurements 
of  A  type  and  B  type  probes  as  well  as 
measurements  of  a  Kiel  probe. 
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Evaluation    of      a      and      <j>      (method) 

The    best   way    to    understand    the    evaluation    of   yaw    and    pitch 
is    to    look    at    it    from    a    topographical    point    of    view.      Each    of 
the    type    A    and    type    B    probes    has    a    unique    calibration    surface 

Cprp    =    Cprp^arp'  Vp^  *    where    arp    and    Vp    are    measured    relative 
to    the    axis    of    symmetry    through    the    sensor    at    the    probe    tip. 

The    calibration    surface    in    this    representation    is    invariant    to 

yaw    and    pitch    of    the    probe    axis    relative    to    the    laboratory 

space.      The    same    calibration    surface    if    represented    in    the    form 

Cp    =    Cp(a,4>)      can    be    written    as 


CP    =    Cp[(ap    +    a       ),     Up    +     4>rp)] 


(5) 


where      a        and      <j>        are    the    probe  tip    axis    angular    settings 

relative    to    the    laboratory    space.  It    is    clear    from    eq .    (5) 

that    Cp    can    be    derived    from      Cp  by    a    constant    translation: 

a      ,       *        on    the      a,*      plane.      As  each      Cp  can    be    viewed 

p  Yp  y      r  r  r  p 


as    a    hill    with    its    peak    at      a        =    0      and      <f> 


0.   The   Cp 


surfaces  are  the  same  hills  with  their  peaks  translated  to 
a  =  ap   and   $  =  $  . 

In  the  present  method  probe  A  is  used  in  three  different 
angular  settings  namely: 


I)  «p  -  o  +p  -  0 
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m)   ap  =  aplll   *p =  ° 

IV)  aP  =  w      *P  =  ° 


and  the  B  probe  is  used  in  a  fixed  position,  namely 


II)   «p  =  0 


'pH 


For  this  case  the  topography  of  the  calibration  surfaces  will 
appear  as  four  hills.   The  three  with  their  peaks  at  points 
(0,0);  (+apIII,  0)  and  (-aDlv»  0)  are  the  translated  hills  Cp   . 
and  the  fourth  with  its  peak  at  (0,  a,,)  is  the  translated 

hill   Cp   n .   Their  contours  of  constant   Cp  are    then  as  shown 
r  r  pB 

in  figure  4. 

Assume  now  that  a  velocity  vector  with  yaw   a   and  pitch 
4)  is  to  be  measured.   These  values  of   a   and  $     will  be 
sensed  uniquely  by  the  probes  at  their  four  angular  settings.   Were 
equation  (4)  single  valued,  the  values  of   a   and   <j,  could  be 
uniquely  evaluated,  as  the  single  intersection  point  between  the 
projections  of  appropriate  lines  of  constant  Cp  on  the  a ,  <f>   plane* 
In  the  present  case  of  the  double  valued  functions,  the  lines  of 
constant  Cp  are  closed  curves  and  more  than  a  single  intersection 
point  do  exist. 


*  Cp  has  generally  a  different  value  for  each  probe  in  each 
of  its  settings.   Thus  the  solution  involves  solving  for  the 
intersection  points  between  projections  of  contours  of  specified 
(but  different)   Cp   values  on  the  different  hills. 
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Fig 


Projections  of    C      =    const,    lines    of    the    four 

calibration  surfaces    (the    center    hill    is    for 

probe   A    at  a      =    0,       <j>      =   0,    the    top    hill    is 

for    probe    B  at      a      =    0,       <|>      =    <£dtt>    the    right 

hill    is    for  probe    A    at      a      =    +aDlIT,       <f>D    =    0 

and    the    left    hill    is    for    probe    A    at      a      =    a    tw> 

^p    =    °- 
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The  situation  is  shown  in  Fig.  5  for  an  example  of  a  number  of 
such  intersection  points.   The  correct  intersection  point,  however, 
is  uniquely  identified  as  the  only  point  through  which  all  four 
Cp  =  const,  curves  pass. 

From  an  examination  of  Fig.  5  it  is  clear  that  the  inter- 
section points  of  three  closed   Cp   curves  projections  only  are 
sufficient  to  identify  a      and  $      uniquely.   However  one  of 
these  three  curves  must  be  that  belonging  to  the  type  B  probe. 

The  details  of  the  numerical  procedure  used  to  obtain  the 
correct  intersection  point  in  the  evaluation  of  a   and   $   , 
are  given  in  the  following  paragraph. 
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Fig.  5:  An  example  showing  a  number  of  intersection 
points.  The  real  one  is  in  the.  second  quad' 
rant . 
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Evaluation  of  the  intersection  point  coordinates  (numerical 
Procedure ) 

The  calibration  surfaces   Cp     are  represented  in  form 

of  a  linear  string  of  values  ordered  into  an  array,  as  shown 

schematically  in  figure  6.   The  numbers  in  Fiq.  6  indicate  the 

position  of  a   Cp    value  in  the  strina.   The  strinq  starts 

r  rp 

with  the  value  of   Cp    at  a  point   (-a   ,  ;  -*„,„•«) 

r  rp  r  rpLo         yrpLo ' 

and    ends    with    the    value    of      Cp    „      at    a    point      (+«      ,,    ;    +  <j>      ..    )  . 

1  rp       r  rpUp    vrpUp 

This  sequence  must  be  kept  and  can  not  be  changed. 

Each  of  the  closed   Cp  =  const,  curves  projections  in  figs. 
2  and  3  is  the  projection  of  the  line  of  intersection  between  a 
plane  parallel  to  the   a ,  <\>      plane  at  a  height  equal  to   Cp   and 
the  calibration  surface.   These  closed  curves  can  be  determined  as 
the  locus  of  the  projections  of  all  the  penetration  points  of  ar- 
bitrary lines  parallel  to  the  a ,  <j,  plane  at  a  height  Cp  and  the  cali 
bration  surface.   Such  penetration  points  are  calculated  in  sub- 
routine "PENPTS",  Fig.  10. 

PENPTS  calculates  the  first  two  penetration  points  of  a 
surface  Z  (in  the  present  case   Z  =  Cp)  by  a  straight  line 
piercing  that  surface.   If  the  surface  is  double  valued  these 
two  points  are  the  only  roots.   The  surface  is  given  as  a  table 
of  numbers  on  a  cartesian  basis   Z  =  Z ( X ,  Y )  (or  in  the  present 
case   Cp  =  Cp(  a  ,  <j>)  . 

The  subroutine  has  the  following  limitations: 

1)  No  roots  can  be  found  on  the  lower   Y  =  const,  boundary 
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Fig.  6:   The  order  of  a  calibration  surface 
array  . 
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2)  The  roots  can  be  found  only  along  a  piercing  line 
which  is  parallel  to  the   X  axis  (but  at  any  height 
above  the   X  ,Y   plane) 

3)  The  table   Z  =  Z(X,Y)  defining  the  surface  must  be 
based  on  a  grid  comprising  lines   X  =  const,  and 

Y  =  const.   The  spacing  of  these  lines  must  not  be 

equal.   In  other  words  the  surface  is  defined  by  a 

rectangular  grid  in  the   X,Y   plane,  from   xml-n   to 

X  a    and  from   Y  ,    to   Ymav. 
max  mi n       max 

4)  Only  a  single  root  can  be  evaluated  in  a  surface 
element  located  above  a  defining  rectangle  on  the 
X,Y  plane. 

5)  Not  more  than  the  first  two  roots  will  be  evaluated  for 
any  pi  e  rci  ng  line. 

6)  The  surface  must  be  monotonic  over  each  rectangle 
(this  is  a  result  of  limitation  4). 

7)  All  X  and  Y  arguments  must  be  given  in  increasing  order 


These  limitations  do  not  restrict  the  application  of  PENPTS 
in  the  present  problem  as  long  as  the  calibration  surfaces  are 
smooth  within  the  element  located  above  a  grid  rectangle. 
However,  the  elemental  grid  rectangles  can  be  reduced  arbi- 
trarily in  size.   If  a  calibration  surface  is  more  than  double 
valued  PENPTS  will  fail.   However  in  this  case  the  probe  yielding 
such  a  calibration  surface  can  not  be  considered  a  useful  instrument 
unless  used  only  in  parts  of  the  domain  where  it  is  double  valued. 
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Z  is  defined  from  X  .   to  X    and  from  Y  .   to  Y 

min     max  min     max 

PENPTS  is  given  the  following  initial  information:   X  values  and 
Y  values  defining  grid  points,  corresponding  Z  values,  the  Y 
location  (YG)  and  the  height  above  the  X,Y  plane  (HT)  of  the 
piercing  line.   With  this  information  PENPTS  searches  for  the  band 
of  rectangles  which  includes  YG  or  of  which  YG  is  the  lower  bound- 
ary (see  fig.  7)  and  then  scans  this  band  from  left  to  right  in 
search  for  penetration  points.   The  scanning  is  based  on  the  geo- 
metry given  in  fig.  7  which  represents  a  particular  element  in 
the  band,  approximated  by  two  plane  triangles. 

Initially  a  coarse  scan  is  carried  out  just  to  detect,  but 
not  to  evaluate,  an  intersection  point.   This  is  done  checking 
for  each  sub-domain  whether  (ZX(I-l)  -  HT) / ( ZX ( I )  -  HT)    0. 
If  this  condition  is  true  a  root  is  detected  and  control  is 
tranferred  to  its  exact  evaluation.   The  value  of  the  root  is 
calculated  after  its  location,  either  in  the  first  (left)  or 
second  (right)  triangle  is  determined  (each  grid  rectangle  is  com- 
posed of  two  triangles).   Equation  6  which  is  based  on  fig.  8 
(for  a  left  triangle)  or  eq .  7  which  is  based  on  fig.  9  (for  a 
right  triangle)  is  used  to  evaluate  the  penetration  point.   These 
equations  express  the  linear  interpolation  of   Cp   in  the  Fortran 
rotation  used  in  this  program. 

XS  =  X(I-l)  +  ABS(ZX( 1-1 )-HT)*(XM-X(I-l ) ) / ( ABS( HT-ZM)  + 
+  ABS(ZX(I-1  )-HT)  )  (6) 
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Fig.  7:   The  geometry  of  a  linearized  calibration 
surface  element  comprising  two  plane  tri- 
angles and  its  intersection  with  a  plane 
normal  to   X,Y   along   Y  =  YG. 
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abs(ht-  zm) 


ABS(HT-ZM) 


Fig.    8:      The    geometry    for    a    penetration    through    a    left 
hand    ( fi  rs t)    triangle. 
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clABS(HT-Z(l)) 


Z(I) 


Fig.  9:   Tne  geometry  for  a  penetration  through  a 
right  hand  (second)  triangle. 
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XS  =  XM  +  (X(I)-XM)*ABS(ZM-HT)/(ABS(ZM-HT)  + 
ABS(HT-Z(I))  ) 


(7) 


Equations  6  and  7  are  invariant  to  the  slope  of  the  cali- 
bration surface  i.e.,  the  slope  of  the  straight  lines  AB  and 
BC. 

Finally  it  should  be  pointed  out  that  when  two  penetration 
points  are  determined  the  values  of  their  abscissas,  X(I),  are 
returned.   When  only  a  single  penetration  point  is  detected  the 
abscissa  X(2)  will  be  returned  with  a  value  of  1000.0.   When 
no  penetration  point  is  determined  both  X(I)  values  returned  will 
have  the  value  of  1000.0.   The  program  logic  is  designed  to 
recognize  these  messages. 

It  was  stated  earlier  that  PENPTS  is  used  to  determine  the 
closed  intersection  curves  projections  on  the  X,Y  plane.   In  fact 
not  the  curves  but  just  the  intersection  points  between  each  two 
of  them  are  required  (see  figure  3). 

To  compute  the  coordinates  of  these  points  the  subroutine 
"INTSCS"  is  used.   It  uses  PENPTS  as  a  subroutine. 

In  INTSCS  a  scanning  procedure  is  carried  out  from  a 
minimal  value  of   Y  (or  $)  YRIN  to  an  upper  value  of  YRUP,  or  a 
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prescribed  10,000  times*  which  ever  comes  earlier.   The  sub- 
routine scans  through  any  two  arbitrarily  chosen  closed  curve 
projections  to  find  their  intersection  points.   In  each  scan 
(I)  up  to  four  penetration  points  can  be  determined,  while  the 
penetration  points  of  the  previous  scan  (J)  are  memorized. 
Together,  eight  penetration  points  can  be  involved.   When  no 
intersection  point  exists  the  geometrical  situation  is  as  shown 
in  figure  11,  while  the  existance  of  an  intersection  point  is 
characterized  in  figure  12.   Subroutine  INTSCS  can  distinguish 
between  the  two  situations.   In  figures  11  and  12  the  case  of 
four  penetration  points  found  in  each  scan  are  shown.   The 
subroutine,  however  will  handle  any  possible  number  of  such  points, 
from  zero  to  four.   A  "no  penetration  points"  is  assigned  an 
abscissa  value  of  1000.0  by  PENPTS  as  explained  earlier. 

When  an  intersection  point  is  detected  its  evaluation  is 
based  on  the  geometry  shown  in  fig.  12.   The  eight  penetration 
points  have  the  following  coordinates: 

AJL(RAJL,YRJ)  BJ L ( RB J L , YRJ ) 

BJR(RBJR,YRJ)  AJ R( RAJR  ,  YRJ ) 

BIL(RBIL,YR)  AIL( RAIL , YR) 

AIR(RAIR,YR)  BIR(RBIR,YR) 


*  The  10,000  scans  are  governed  by  a  program  constant  in  the 
line  "DO  140  ...."  and  can  be  arbitrarily  varied. 
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Fig.  11:   The  geometry  for  two  successive  scans 
when  no  intersection  point  exists. 
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Fig.  12:   The  geometry  for  two  successive  scans 
when  an  intersection  point  does  exist 
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The  left  intersection  point  illustrated  in  fin.  12  is  the 
intersection  point  of  the  lines  aa  and  bb,  each  described  by 
its  equation: 

for  aa   Y  =  (AA)X  +  BA  (8) 

for  bb   Y  =  (AB)X  +  BB  (9) 

The  constants  AA ,  BA  ,  AB  and  BB  are  given  in  equations  10  to  13. 
AA  -  (YRJ-YR)/(RAJL  -  RAIL)  (10) 

BA.  =  YR  -  AA*RAIL  (11  ) 

AB  -  (YRJ-YR)/(RBJL-RBIL)  (12) 

BB  -  YR  -  AB*RBIL  (13) 

The  coordinates  of  the  intersection  point  to  be  calculated  are 

X  =  (BB-BA)/(AA-AB)  (14) 

Y  =  X*AA  +  BA  (15) 

These  relations  are  true  for  a  left  hand  intersection  point. 
Analogous  equations  are    true  for  a  right  hand  intersection 
point.    In  this  algorithm  the  straight  lines  aa    and  bb  approxi- 
mate the  curved  lines  connecting  AJL  with  AIL  and  BJL  with  BIL 
or  similar  lines  on  the  right  hand  side  of  figure  12.   The  error 
introduced  through  this  approximation  is  reduced  as   DAR  =   A  $ 
is  reduced. 

It  is  possible  that  an  intersection  point  is  identical  with 
AIL  and  BIL  or  AIR  and  BIR.   This  case  is  defined  as  "direct  hit" 
The  program  is  designed  to  detect  such  a  direct  hit  and  evaluate 
the  corresponding  intersection  point  accordingly. 

The  above  algorithm  works  perfectly  as  long  as  the  two 
closed   Cp  =  const,  curves  are  far  from  being  tangent.   But  in 
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practice  a  situation  of  almost  tangent  curves  can  arise  when  X 
(or  a)  is  very    small.   In  this  case  the  preceeding  algorithm 
will  fail  and  must  be  replaced.   The  geometry  of  this  situa- 
tion is  described  in  figure  13.   This  situation  is  identified 
by  INTSCS  and  the  intersection  points  are  then  evaluated  assuming 
that  they  are  intersections  of  two  circular  arcs.   When  the 
curves   Cp  =  const,  are  almost  circular  this  approximation  does 
not  lead  to  unacceptable  errors. 

It  was  stated  earlier  that  INTSCS  scans  from  a  minimum 
value  of   <j>   to  a  maximum  value  of   $  with  prescribed  steps 
£<}>,   as  shown  in  figs.  11  and  12.   This  direction  of  the  scan- 
ning is  used  when  the  intersection  points  of  the  curves   Cpr 
and   Cpyr   of  fig.  5  are  evaluated. 

However,  in  course  of  the  reduction  of  the  measured  data, 
scans  in  the  direction  of  a      in  steps  of  A  a   are  also  neces- 
sary to  evaluate  the  intersections  of  the  curves   Cpr   and   CPttt 
or   CpT  and  CpTy.   INTSCS  is  designed  to  carry  out  this  task 
as  well.   To  do  this  the  calling  statement  for  INTSCS  is  appro- 
priately changed  as  will  be  explained  in  the  next  section.   To 
be  general  enough  INTSCS  is  not  written  in  terms  of   a   and   $  or 
X  and   Y   but  rather  in  terms  of  general  arguments.   The  best  way 
to  understand  INTSCS  is  to  compare  its  general  arguments  to  physi- 
cal quantities  by  means  of  the  calling  statements.   In  fig.  14  the 
flow  diagram  in  INTSCS  is  given  in  terms  of  the  general  arguments. 

INTSCS  returns  to  the  main  program  (fig.  15)  the  (a,<}>) 
coordinates  of  two  intersection  points  between  the   Cp  =  const, 
curves  specified  in  the  calling  statement.   Let  us  now  follow 
the  way  in  which  the  main  program  is  designed  to  utilize  INTSCS 
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Fig.  13:   The  geometry  when  two  closed 
curves  are  almost  tangent. 
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for    the    evaluation    of      a      and       <t>     of    the    velocity    vector,    as 
well    as    static    and    total    pressures. 
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Evaluation  of  the  Velocity  Vector  and  Pressures  from  the  Probe 
Signal s 


Fortran  program  VDR  was  written  to  evaluate  the  velocity 
vector  from  probe  measurements  of  pressures.   The  program  is 
shown  in  Fig.  15. 

At  line  1410  INTSCS  is  called  to  scan  curves  I  and  II 
for  possible  intersection  points.   Scanning  can  be  carried  out 
in  the  direction  of  the  ordinate  only,  with  the  calibration 
curve  matrices  compiled  exactly  as  shown  in  fig.  6.   This  limita- 
tion is  imposed  by  the  way  PENPT5  is  constructed.   In  the  case 
of  the  intersection  points  between  I  and  II   ZA  and  Z  B  are 
scanned  without  difficulties  in  the  direction  of  the  ordinate 
which  is  <j>   as  shown  in  fig.  12  and  returns  with  the  coor- 
dinates of  the  two  first  intersection  points,  points  1  and  2 
of  Fig.  16.   They  are  ALF1,  PHI!  and  ALF2,  PH 1 2 .   In  line 
1500  INTSCS  is  called  again  to  scan  curves  I  and  III.   Now 
scanning  has  to  be  carried  out  in  the  direction  of  the  abscissa, 
a  task  for  which  INTSCS  was  not  designed.   To  overcome  this 
problem  the  calibration  curve  matrices  are  used  in  a  transformed 
form  such  that  the  previous  abscissas  are  now  ordinates,  or- 
dinates  are  abscissas  and  the  internal  structures  of  the  Cp 
arrays  are  accordingly  modified.   This  transformation  is  carried 
in  VDR  in  the  section  between  lines  320  and  450.   Comparison 
of  the  calling  line  1500  to  the  previous  calling  line 
1410  shows  very    clearly  how  the  various  arrays:  original  and  trans 
formed,  are  used.   The  coordinates  of  points  3  and  4  of  fig.  16 


43 


IRAKI  Mil 

CALIBRATION 

MATRIX 

Sfl  COMSTAMTS 

A«0  FLACS 


READ  PI 


5 


INPUT 
NX,  NY,  XA,  XB,  YA,  TB,  ZA,  ZB 


©- 


©- 
©- 


©- 


YES 

<WSIM-J> 

^Tno 

READ 
3HULATEO 

EXPERIMENT 
DATA 

<WSIM  -p> 

YES 

READ 
PROSE 
SIGNALS 

' 

r 

READ  PSIM 
C INITIAL.  GUESS 

OF  STA 

C  PR) 

ITT  - 1 .  1COPS  ■ 

ISCAN  -  ISCAN  ■ 
PSIN-PSIN>5  0 


PS  •  PSII 
-P  PSC  •  PSIN 


CORRECT  PS 

IF  CP  VALUES 

FICfEO  MATRIX 

RAMCE 


■j.  a  milt  lOOO 
TINES- 


READ 
ADDITIONAL 
DATA  AID 
6ECIN  MAIM 


JTES 

IT  >10 

"nc 


TIT  - 1  IT  +  1 


CALCHATE 

NEW  PS 

SET 

PSC  t  PST 


-0 


-0 


0 


EVALUATE 
IEPS1,  ALFI,  PHII, 
•IF?,  PHI2.IEPS? 
ALES.  ALF4.  PHI3, 
PHI4 


SELECT 
RICHT  a   I  * 
F»LF,  FPHI 


CALCULATE 

NEI 

STATIC  PRES 

PS! 


EVALUATE 
EPSPSC  FOR  CASE 
Of  ALMOST  TAN- 
CENT  CP  CURVES 


PRINT 
FALF.  FPHI. 

PS«,  PT.   AM 


TO  PASS 

THIS  POUT 

CT  MUST  HAVE  SEEN 

REMOVED  FROM  COL  I 


PRINT 
VIR1IM,  CPC 


IEAO 

ADDITIONAL 
DATA 
KCIN   ACAIN 


t  I  ijUP.E. 


FLOW    DIAGRAM    OF    vDR 


* 


—  OL 


Fig.  16:   Selection  of  the  proper  intersection 
point  in  VDR. 


45 


are    now    returned    to    the    main    program    of    VDR,    they    are    A  L  F  3 ,  P  H I  3 
and    ALF4,PHI4. 

We    are    seeking    a    single    intersection    point;    the    one    repre- 
senting  yaw    and    pitch    as    sensed    by    probe    A    in    positions    I    and    III 
and    probe    B    in    position    II.       But,    as    the    calibration    surfaces 
are    double    valued   we    are    now,    unfortunately,    in    the    possession    of 
four    points.       The    solution,    however,    is    physically    unique.       Only 
a    single    velocity    vector   exists    in    reality    and    its    yaw    and    pitch 
are    included    in    the    four    intersection    points    evaluated.      Were 
both    the    measurements    and    the    numerical    procedure    absolutely 
accurate,    two    of    the    four    points    would    have    been    identical.      But 
this    is    not    the    case    in    reality,    instead    of    a    sinole    point,    two 
points    close    to    each    other   will    be    detected.      Therefore    the 
average    coordinates    of    the    two    of    the    four    intersection    points 
returned    to    VDR    from    INTSCS    which    are    closest    to    each    other 
are    selected    as    the    measured   yaw    and    pitch    angles.       In    the    ex- 
ample   shown    in    fig.    16    this    will    be    the    point    dividing    the    dis- 
tance   between    points    2    and    3. 

The    calculation    is    now    at    the    point    at    which     ox       and     <j.       are 

temporarily    known    (see    fig.     3).       Using    a    and    4>    a    new 

CpIV    (for    probe    A    in    position    IV)    is    computed    by    linear    interpo1< 

tion    using    subroutine    INTPLT.      With    this    new    value    of      Cp      TW      and 

Krpl  V 

with  PPtw  a  new  P   is  computed.   If  this  new  value  of  P   is  close 
enough  as  determined  by  EPSPSG  to  the  guessed  value,  or  to  the 
value  of  P   in  the  previous  iteration,  the  data  reduction  for  the 
particular  point  in  question  is  terminated. 
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The  relative  difference  between  present  and  previous  Pr 
is  compared  to  the  convergence  criterion  EPSPSG.   This  criterion 
is  evaluated  by  an  empirical  function  determined  to  give  best 
compromise  between  accuracy  of  results  and  computer  time  required 
until  convergence  is  achieved.   When  the  calculation  follows 
normal  routine  the  convergence  criterion  is  given  in  line 
1710  as  function  of  PT  .   When  the  routine  for  almost  tangent 
curves  is  used  during  data  reduction  either  or.  the  right  or  the 
left  side  a  different  empirical  function  (lines  1720  or  1730) 
i  s  used . 

If  the  convergence  criterion  is  satisfied  results  are 
printed  out  and  data  for  a  new  measurement  point  is  read  for 
reduction.   If,  however  convergence  is  not  reached  a  new  static 
pressure  for  a  next  iteration  is  evaluated  (line  1780). 
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Convergence  and  accuracy 

The  convergence  of  the  present  iterative  procedure  is  not 
ensured  in  all  possible  cases  of  data  sets.   The  program,  however, 
is  adjusted  to  converge  in  most  of  the  cases.   Similarly  to  other 
iterative  computation  method  some  experience  is  required  to  achieve 
convergence  when  the  calculation  does  not  converge.   In  this  para- 
graph the  principal  factors  affecting  convergence,  computation 
time  and  accuracy,  which  are  abviously  coupled,  are  pointed  out. 

The  convergence  and  accuracy  of  reduction  of  a  data  set:   Pj, 
PPA,  PPB,  PPC,  PPD,  PSIN  (or  analogous  set  in  the  experiment 
simulation  mode)  depends  on  the  following  factors,  which  can  be 
varied  by  the  user. 


1)  Coarseness  of  calibration  arrays  ZA,  ZB  and  their 
linearity.  With  coarser  arrays  convergence  problem 
will  increase  and  accuracy  of  results  decrease. 

2)  Probe  settings   YPB,  XPC,  XPD  (with  the  following  probe 
settings  being  fixed  and  not  variable  YPA  =  YPC  = 

YPD  =  0.0;   XPA  =  XPB  =  0.0).   The  optimum  setting  is  about 
^20°  to  25°.  Too  small  values  reduce  accuracy,  too  big 
values  cause  convergence  problems  and  probe  tip  flow 
separation. 

3)  RELXPS,  the  static  pressure  relaxation  factor.   The 
smaller  this  factor  the  safer  will  convergence  be 
achieved.   Computation  time,  however,  will  increase. 


4) 


The  constant  5.0  in  line  1050  of  the  program.   The 
dimension  of  this  constant  is  kg/m  .   If  too  big 
results  can  be  lost  and  if  too  small  computer  time 
will  be  growing.   The  variation  of  this  constant  should 
be  coupled  with  an  appropriate  modification  of  the 
constant  in  line  1010  of  VDR  (item  8  in  this  list). 
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5)  The  constants  in  the  evaluation  of  EPSPSG  (lines  1710 
to  1730).   The  smaller  EPSPSG  the  more  accuarte  the 
results  in  expense  of  increased  computation  time  and 
reduced  convergence  safety. 

6)  The  constant  scanning  step  DAR.   The  smaller  DAR  the 
more  accurate  the  results,  but  too  small  values  can 
cause  complete  loss  of  results.   Computation  time  in- 
creases with  reduced  DAR. 

7)  The  constant  10  in  line  1360.   This  constant  governs 
the  number  of  P   corrections. 

8)  The  constant  1000  in  line  1010.   This  constant  governs 
the  number  of  s  cans  . 


If  convergence  is  not  achieved  in  a  particular  case  varia- 
tion of  each  or  of  a  combination  of  the  above  values  will  always 
enable  convergence. 
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Evaluation  of  the  Computation  Time 

By  deleting  the  letters   CT  from  column  1  and  2  of  lines 
20  to  70,  1930  to  1960,  1980,  2000,  2010  the  actual  computation 
time  as  well  as  CPU  time  are    evaluated  and  printed  out.   The 
following  statement  prior  to  execution  is  required  in  this  case:* 

GLOBAL  T  SYSLIB  SSPLIB 

This  option  is  useful  for  adjusting  the  constants  affecting 
convergence  as  to  optimal  compromise  between  accuracy,  ease  of 
convergence  and  calculation  costs. 


*  When  the  program  is  run  on  Naval  Postgraduate  School  IBM  360 
sys  tern. 
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I  nput 


1)  Calibration  arrays  have  to  be  input  in  the  following  manner: 
first  values  of   NX   and  NNY  ,  second  values  of   X   in  rising 
order,  third  values  of   Y   in  rising  order  and  fourth  values 
of   Z   in  the  order  shown  in  Fig.  6.   Two  calibration  arrays 
are  read:   first  the  array  of  type  A  probe  and  second  the 
array  of  type  B  probe. 

2)  If  the  program  is  run  in  experiment  simulation  mode   PT 
and   PS   are  read  in  Kg/m2  and  subsequently   o   and   $  . 

3)  If  the  program  is  run  in  data  reduction  mode  the  measured 
values  of   PT,  PPA,  PPB,  PPC   and   PPD   are  read  and  then 
the  guess  of  static  pressure   PSIN.   All  are  read  in  Kg/m2. 
PS  I N   has  to  be  lower  than  the  actually  existing  value  to 
ensure  convergence  of  the  calculation.   Too  low  a  value 
will  cause  waste  of  computer  time. 

Calibration  arrays  are  read  from  a  disk  space  on  which  they 
are  stored.   The  following  statement  prior  to  execution  is  re- 
quired for  successful  reading: 

FILEDEF  02  DSK  NAME  XX 
here  NAME  is  the  name  of  the  file  on  which  the  calibration  data 
of  probes   A   and  B   is  stored  in  proper  order  and  format,   XX 
i  s  a  two  digit  number . 

Simulation  or  measurements  reduction  input  is  read  in  the 
normal  way  using  the  terminal  keyboard  or  punched  cards. 
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Output 

a,  4),  P  ,  Pj   and   Ma  are  printed  out  according  to  lines 

1820  to  1920.   This  output,  however  is  not  sufficient  when  the 

logic  of  the  computation  is  to  be  followed  either  to  examine  the 

execution  of  the  program  or  for  debugging.   Two  programs  with 

additional  output  are  given  in  appendices  2  and  3  and  can  be  used 
for  this  purpose.   The  first,   SWVDR   gives  short  additional 
output ,  namel y : 

1)  When  PSIN  is  successively  increased  automatically 
by  the  program  to  ensure  convergence,  the  values  of 
PSIN   are  printed  out. 

2)  IIT,  FALF,  FPHI   and   PSN   are  printed  out  at  the 
end  of  each  iteration  prior  to  the  convergence  test. 

The  second  program  WVDR  prints  more  detailed  information.   All 

additional   WRITE   statements  in  this  program  are  numbered  with 

three    digit    numbers    starting   with    9. 
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Concl usi  ons 

The  following  conclusions  are  drawn  on  the  basis  of  experi- 
ence gained  in  running  the  program  with  various  data  sets  and 
various  types  of  calibration  surfaces. 


1)  Optimal  probe  settings  for  ease  of  convergence  and 
high  accuracy  are:   YPB  =  25°,  XPC  =  25°,  XPD  =  -25°. 
It  is  therefore  suggested  that  a   B   type  probe  with 
25°  pitch  will  be  used,  and  that  the  A  type  probe 

be  rotated  to  +  25°. 

2)  Convergence  and  accuracy,  as  well  as  computer  time 
efficiency  are  improved  when  the  calibration  surfaces 
of  both  probes  are  not  flat  at  their  peaks  but  are 
rather  rounded.   It  is  therefore  suggested  that  a  new 
probe  tip  geometry  be  considered.   A  spherical  tip  with 
a  central  pressure  tap  is  recommended.   To  prevent 
damage  to  the  sensitive  transducer  located  behind  the 
pressure  tap,  and  in  order  to  improve  the  frequency 
response  of  the  probes  it  is  suggested  that  the  volume 
between  the  pressure  tap  face  and  the  transducer 

be  filled  with  an  appropriate  liquid  that  will  not 
affect  the  transducer  negatively.   In  this  case  the 
opening  of  the  pressure  tap  has  to  be  sealed  with  a 
very  thin  low  inertia  membrane. 

3)  It  is  probably  possible  to  modify  the  iterative  pro- 
cedure such  that  safe  convergence  can  be  achieved  also 
when  using  the  scheme  of  Fig.  2.   If  this  can  be 
achieved  the  Kiel  probe  will  not  be  necessary.   An 
effort  in  this  direction  is  suggested. 
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APPENDIX  I 
LISTING  OF  VDR 
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U  I  "IN 
X,  NY 


i.  j  )  j    \  M  1  i 

1  =  1   ,f\Y 

1 3  c  j    r a (  I J 

1=1 t    l< 

J  =  l  ,  u 

( 1 1 

NX 

x~m 

IY 

Y  i(  1  J 
'IX 
\Y 


13C) 

1=1, 


RtAJ(?  ,119) 

no    121   1  =  1, 'ix 

121  REACU. 
UJ    122 

122  ^65,0(2, 
CO    123 
0  J    123 

123  R£^C(2,1^JJ 
CU    22i     1=1 

221  REaD(2, 
CQ    22^ 

222  FEAG(2,13C) 
C3  223  1=1, 
CO    223     J=i , 

223  READ(2il2CJ 
C    TRANSFCR^    :  ALIRRATUN 

OL    224     1  =1  ,.'W 

224  XAX( 1)=  i4(  i  ) 
OG    2^b    1=1, 'IX 

225  YmM  i  J  =xa(  i) 
03    2  2  c-    I  =  i  ,  M  f 
CC    2Z6    J=i , N< 

22t  ZA>{  I  ,J)  =ZA(  j  ,1  ) 

03    ZZ7     1=1, 'IY 

227  XdX(lJ=YH(!) 
03    22t    1  =  1,  IX 

228  YE!X<  I)  =  XE(  1  J 
DC    229    1=1, NY 
03     2  <9    J=l,'lx 

229  Z6X(  I,  J)  =  Zo(.J,  I  ) 
C     FRCEl     icIfiN^G 

YPE=55  . 
>P  A=G.  J 
XPC=20. 
XPL=-2J. 
C    PR3GRA*"    C  C  s\  5  I A  M  I' S 
YRLH  =  -jJ  . 
»KLF=3  J. 
YR  1M  =  -3C. 
XRLC=-jJ. 
>RIP=3G. 
>R  1N=-jC. 
ICCPS=1 
NOCOP?=l 
NO S I*=l 
UT  =  1 

iscai\  =  i 

9CLXPS=J  .5 
C    EXPE-il  Nc  M     Si  ILLATION1 
<«0  1         F3RMATIF1S.4) 

193  Fl^C  (i,  =tJl  )     JT 

fcRlTE(6,4Cl)     PT 
1FIMCS1 -1.hj.t  J    3i 
RtAO(5,40l)     PS 


D'->7\ 


J  i 


J) 


•'A  TRICES 


73     500 


IN  roju  la 

INTG0020 

INTlCCIO 

I  n  ojj4j 

INTO  5  0  '30 

INTO  JuOU 

INT 0  30 7 J 

INTCJCeO 

IN  'C^OO 

INfOJlOO 

INTCC1  10 

INTO^lzO 

IN  1001 j  ) 

IN  TOG  140 

IN  T0U15O 

I N  T  C  G I  c  G 

INT0J1  7'J 

I N TOJlao 

INTlCISC 

INT0J2 J  J 

1NT0U210 

IMYCQ220 

INT0JZ30 

1NTCG2  4G 

INT002^0 

INTCJ260 

I  U0J2  70 

INTOGZdO 

INTCC2S0 

INTCG2C0 

I  NT  00 3  10 

INTCC320 

INT0O33O 

INTC0340 

HTCG35Q 

INT0J36J 

INTCC3  7C 

INTO  J  3^0 

INT003  90 

I\|TCG4C0 

INT-JJ4  1J 

INT0G42J 

INTCC433 

INT0044J 

INTG'J4t.J 

lNTCC4cG 

INT0J47J 

INT JO 430 

iNTOO'iyO 

iNTJOiJOJ 

INTCC510 

INTC052J 

INT00b3J 

INTLC54C 

INT0O350 

INTuo5oO 

INTCGi70 

INT005  6J 

I N  T  C  C  5  S  C 

IN  TOD  6 JO 

I  NToOCiO 

IMTCC62J 

IN  TCC6 JO 

INT006  43 

INTCC65C 

INTJ0t6J 

INTOJ6  70 

iNTCvCto 

INTO J69J 

INTCC7C0 
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FILE:    VJt- 


FHRTPAM       PI 


N 


<'iim- 


k>RI  TE(  c.Oll     P", 

READ(5,VJ1)     iLf 

W  8 1  T  £  l  c.  ,  t  j  I  )     ALFA 

REAO(  5,  ^  C  i  J     °rii  I 

wR  ITE  <o,  +Jl  J    Whii 

CALL     1  <  TPl  T(  aA  ,Ya  ,ZA  ,ALFA  ,  Prill  »NX, NY, CP£) 

alk  =  vf\-ap; 

CALL     I^TSLT  (x  afy  A,ZAf  ALFC»  PHII  ,NX,  *JYf  C«>C) 
AlFO-^  lFa-aPJ 

tZA.ALFOf  PHII  ,MX,  JY.CPJ) 


CALL     J  ITPLT  UA,y  £ 
Pril  a=FnI  I  -YPd 
CALL     IN|PLTU3,YF 
FOfN=PT-Pj 
PPA-CPA*FJ YN+PS 
PP3=CP  -?*PJY^+P5 
FPC=CFC*FJYN*PS 
PPO  =  CPJ*PJY\«-Pb 
C    PEAO    McASUBt'-IEMTS    CAT 
I  F  i  N  c  S  i  ■' .  t  J  .  1  j     C  C 

WRIT  £  (o,  ijo  J    ?->A 
At  AJ  15  ,135  J     PP3 
WR  1TC<  c,  Ufc  )     P»  * 

Hc/C (5 ,1 iO I       PPC 

k»RlTE(  L,lJ'wJ     PFC 
3E^C(5, li5)    •>  °  D 
WRITEU,  l.io)     PPJ 

5ol  IE  a0(  5,  13; )     PS  1  4 

WRITE  (o,  ijoi     P^  IN 
IF(iSC^N.tw.l)    o^ 

18C  IF  (  ISCV4.0  I.  1JJJ) 

IlT  =  i 
ICUPS^l 
ISCAN=1SC  VJ+-1 
FSi.\=PSifn-5.0 


,  Z  b  ,  A  LF  A  ,  PH  I  el  ,  NX  ,  NY  ,  C  P  R  ) 


TC    501 


TO    131 

J  J    TO    30  5 


181  PS=°3  1m 

F  SC=  F£l  n 

psr=ps 

C    CALCLJLATt^     >RcSSU3 


30C         P0YN=PT-F5 

Cpa=  (  pp4-p  j  j/p:)y;>! 

CFb  =  (  FFd-PS  J/PJYN 

cpc  =  ipp;-p j)  /pdyn 
:fl>  (fpc-psj/^cyn 


COEFFIC  IENTS 


COFFECTS  Pi  ;S^L-lPTiCN 
IFfMflPi.cJJi  ,J 
2A*'AX=-100u0  . 
CAlL  M  *Xa->Y  (  ZA  ,Z  a*A 
IFICPA  .J  f  .*.*  <UX  )  Pi 
IF  IlFC.  J.i^uA  ) 
IF  {  CPO.o  I.  Z  A  MAX) 
2E^flA=-  lJ^J  . 
CALL  KAXA-Y  4^3  ,Z 
IMC™. 31  .Z  3  1AXJ 
Za*U'=  10  JJO  . 
CALL    M  NA  i-V  ii.ntl^ 

ificpt  .I.T.Z4--.  rii 


Tr 
1  ] 


ENSURE 
3U1 


CF    ARE     IN    CALidRATlJN    RANJE 


X  ,NX  ,NY  ) 

_C  =  1.  Jl-^S 

CSC=1  .01  *PS 

PiOl.dl* PS 

[v^x   ,NX  ,  NY  I 

ps; =1.01* ps 


i 


IF  (CFC. Lf . i j 
IF  (CPO.^T.  ZA 
ZdM"  =  1JJJU. 
CALL     Mi|,irVUr 
IFICPB.L  I.Z"1.  <I 
IFIPSC.ES.PSI  i 
IC0P5=1  CCPS-H 
IF(  lCr,J0  .JT.10J 
FS^PSC 

pst  =  p:>c 

CO    TO    jJJ 
CALCLLA  lhS     T/J    AcF 


MY  J 
PoC=  C.'jv*0  S 

p->c=o.?y*p$ 

F3C  =  C.  vc/-PS 


,NX 


iZE 

1  j       ! 


!  ,NX  ,i\Y  ) 

=  C.   V9*P 

JJl 


GJ    TO     idO 


t     Fril     ANJLeS    EOF    'i 


1  I 


INT  J 
I  MO 
IN  TC 
IN1C 
I N  T  j 
INTO 
i  NT  0 
INTC 
IN  TC 
INTj 
I  i  TC 
INTJ 
INTJ 
INTC 
INTJ 
INTC 
IN  TO 
INTJ 
IN  TC 
INTJ 
INTJ 
INTJ 
I  NT  0 
INTC 
INTO 
INTO 
INTC 
INTJ 
INTJ 

in  r  g 

INTJ 
i  vjTO 

INI' J 

INTJ 
INTC 

INTJ 
INTC 
INTJ 
IN  TJ 
INTC 
INTJ 
INTC 
INTJ 
INTJ 
IN  TC 
INTO 
INTO 
IN  TO 
INTJ 
INTC 
INTC 
INTu 
IN  TC 
INTJ 
INTO 
INTC 
INTO 
INTC 
INTJ 
INTC 
IN  TC 
INTO 
INTC 
INTJ 
I  NT  0 
IN  TC 
INTC 
INTu 
INTJ 


J  7  1J 

CI  JJ 
j7AJ 
J/5J 

0  I uJ 
J/  /J 
J7rtJ 

u  7  *•;  j 

JdOJ 
^fi  1J 
J  HZ  J 

J.IJJ 

JJ*>J 

ICUu 

H70 

J  J  Ju 
C-SJ 

■;V  JJ 

0910 
J 9  20 
OS  J  J 
CSAJ 
u9b0 
J<?oJ 
C970 
J'i  J  J 
J94J 
ICCJ 

u  lJ 
1C?J 
10  30 
1040 

iceo 

10  7J 
1CHJ 

1  J  90 

I  Uu 

11  U 
11Z0 

II  JJ 
11  40 
1150 

I  1 1: 0 
U7J 

II  JJ 
1  100 
1Z0J 

It  u 

1^20 
itJJ 
l^C 
1/30 
I^jJ 

i :-  70 

1<_0J 

12C.C 
13J0 
i  J  iJ 
1220 
UJJ 
1340 
Ij^J 
1  JoJ 

I J  JJ 

1  JEJ 
1  jOJ 
14UJ 
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FILE 


VUR 


FCRTKAN       PI 


3C1  CALL     I  U33S(  YPA  ,YPl,  YULC,  Y  FUP,Y'riv!  , 

1CFA  ,CPef  ALh  ,  JHi  ,  a4,  yA.Z*,  , 
2Xt),Yfc,te,\fiX,Y^A,^/lX,XtX,YtJX,ZeXf-)X,r!Y,IFPS) 

1 1£  P  S  1=  i  2  .^  S 

*LH-<LHl  J 
ALl  ?=/Si_»-  (2) 
P  H  1 1  =  P  rl  I  (  i  i 
PHl<>  =  FnU2  J 
C    CAcCLuAKb     IaJ    «lF     R     FU     ANGLES     F  D  K    M' 


1  11 


C    Sfcl 


lCPAtCFC,f.iI,ALJ-,X/X,Y.iX,z:^X, 
2>A>,  YA  a  ,L\  /,  A4  ,  ITm  ,£A  , /A  ,  Y~  ,ZA  ,  NY  ,' 

IEPS2=  It5", 

ALf  2=alF  (  1) 

ALF4=ALF{ *) 

FHli^f.il  m 

PHI 4=FHi ( 21 
LECTS    THl     P-MPtR     Al  F     R    Phi     \i 

FALFMt-LFl  »ALr2  ttcFJ  +  ALF- 
1UF1  )-i  ^U  i.  I  Ai_F  I  ,  A^F2  ,-.  lF.j 

FPF  1=  (  (*>h  11  +  >F  i_  +  Jh  I 
1-A'U,U(  Pl-tll  ,PHi2  ,P,iI 
C    CALCULATES      I  Ed     r>i    F*CV 

ALFCf^=F  aLF-X  PJ 

CALL     I  l>  I  PL  I  (  Xm,  if 

4LFCN=FiLF-X  JC 

CALL     IMfLTl^,y4,/^,;LFC'l,FPh  it  MX,  NY,  CPU 

PSN  =  (PPC  «*.  ^l)-PI'D*C^C  l/(  CPJ-C  PC  1 
C  CCNVEKLtNLc  Tl^TS 

EP^PS^i  il  I  Fi^-FSl / P  S J 

EPoPSC=  .'J  JuCjoj  J  (J  C  V  9  <  (  PT*  -2.  C1-.  000  D 00 98* PT 

IFUEFS1.lu.2J     tPSPS  ^  =.000  0000  J009*(  P7  **2  .J  )-.  0000  00  c^T 

1F(IFP$2.EQ.2)     EPS.3  SG  =  .  o  COCUOOc  Gu9*  (  PT**  2.  C)  -.  0  0  00  0  0  8*  DT 

IF  (fcPSPS  .  JT  .EPSPoG)     LC    TU    310 


NX  ,  IEPC.  ) 


iOLE3    HOT    IDF     THE    FJUK     VALLEY 
j- V4AX1  (ALF1,  ALF2,  ALF  j, 
.  Li-  4)  )  /2 

♦  OhI-t)-VMM.(r,HIltPHI^»?'Uj,PHl4) 
t°HL'*n  / 


OAT  A    OF     P  13  IT  ION     IV 
,1m  ,  ALFD  -J  ,FPril  ,  NX  ,  ,\Y»C  PO) 


310 


303 
171 

172 

173 

174 


175 
CT 

CT 

CT 

CT 

450 

CT 

451 

CT 

CT 

3C5 
30  7 


GO    TC     3C3 
1 1T=  1  IT  +  1 

IF  (  I  IT. 01. LO)     jl    TO     130 

P:>  =  PS**ELAPi*I  Pjh-PS) 

FSC=PS 

PST=FS 

GO    TCI    jOu 

UkITL(of17i)     FALF 

FJk  4AT(  •  o'  ,  •     VtL'.CI  TY     YAW 

^RlTE  <6  ,1/ ^)     FFn  i 

FOF.."ATl  ■     VlLCCI  TY     PI  TCn     IS 

^riL(o,l/;J     '3^ 


FCFN^Tji  ,T;Tk 
ViRITE  (  o  ,  1  7  hJ  pr 
FCPVAT(  '     TOTAL     P 


F-cS"UaC     IS 


IS'     ,F15.2) 
'      ,F15. 21 
,  F15 .3  1 


E.oL'E     IS'     ,F  15.j) 
AP=StFT  ( {  I  FT /PS. \  ;**0  .26^71 4-1  .0  )  *5 .0  1 

wrkeu,i/3)    v-i 

FOF'^Tl*      -1  »CF     ••10SRER     IS'     ♦  F  15.  5) 
CALL     I  ALL  ?,\(A) 
VIRT  !»*=  {  M5  )-  iA:>~  )/  ittiuj 
CPL-lA(tij-iAf.  l)/  7ft3  0  0 
rf*ITE(c,4;i^)     V  1  :,  1  i  1 
FO^/.T(«     vinTLUi.    tJm£    13'     ,F15.b) 
1»R1  TE  (5  ,45  iJ     CPU 
F3RMAT{  •     T  )T\L    CPU    IS  •     ,F  lb.  5) 
I  A5G=Al3 1 
IA  6C=m(  61 
GO    T"i    19  2 
«.U  Tl(  'j  ,30  7) 

FCPMATt  «     jCAri^C     1000    TIFFS') 
oO    TC    19  J 
EHi) 

SUE  ROUT  iNE     IVTiCit  NRP1,  \R02,A?.LT,  \P-UP,  A?  IN  , 
lHTl  ,HT2  ,iol  ,  Ri3^  ,X1  ,Y1  ,/l 


2X2 


,Y 

CI^ 


liit  <:  »  a  j  ,  f3  ti.  .;  ,  a  *t  ,  Y  '♦ ,  t '»  » K  1 

E  K  5  I  1  J    J.  E  3  i  (  iJ  I  ,  ■-,  E  3  t  (  i  J  1 ,  A 


'4  ! 


I  EPS) 


Eo  I  1  J)  ,  XI  (4j)  ,  r  1(  40  )  , 


INT  0  li  10 
IMTCltJO 
IN  TC  142C 
IN  TOit^O 
INI  01^ DO 
INTO  l-'«fcO 
INTO  It  70 
1NTC14EU 
IMC14S0 
INTO  1500 
1NTC1510 
IMT0  1520 
IN TQ  1530 
I:lTCl340 
INTO  15 50 
IN  TLLbcC 
INT  01  3  I  ) 
IMT015HJ 
INTC15S0 
INTUIT  JO 
1  h  7  C  1  6  U 
INTO  11  2 J 
INT  01 63  0 
INTCU40 
1NTC1c5G 
INTCUtO 
I N  T  C  I  <j  7  C 
INV016dO 
I  Al  T  C  1  c  5  0 
INTL1  7CC 
I  N  T  C 1  7  1 0 
IN  TCI  720 
INTOl  7j.D 
INTC174J 
INTC1750 
INTO  17 60 
1 N  T  c  L  7  7  C 
INT  01 760 
IN  TO  17  90 
INT  C  La  CO 
INTO  id  10 
IN  TCI  E20 
INTOliiJO 
1NTCIE40 
INTCii5u 
INTO  18 oO 
IMTC1E7C 
INT  01  d60 
INT  Gib  90 
IN  TCiSCO 
INTO  IS  U 
I  JT0192O 
IN "l  CI  9  j  J 
INTC1S,4J 
iNTClStu 
INTJi9oJ 
INT  Clc.  7C 
INT01980 
INTC1c9j 
INTC2CC0 
INTL2Clo 
1NT02O20 
INTO2030 
It\lT02u40 
INTC2C50 
INTJ2060 
IN  TC20  13 
IN TC20JU 
INT  020  9  J 
I NTO^l  CO 
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Flit:    VL>- 


FURTRA*       31 


N 


1Z  1(  4C.43)  ,X2(  40)  ,  iH'iOt 
22->(sO»-«0),Ai(40j,r4{*C),Z-*litO»-iv/) 

RtSUlJ  =  i  MJ. 

Rt  5  1  (  i)  =  >jJoO. 

FE"^2  UJsWJJ. 

RE  StiZt  =S0  10. 

IEPS=1 

IS  t  =  l 

ISP  -  1 

CAR=1. 

4R2-<-A*ir* 

go    r:   15  ^ 

150  <P2r-<-/-<2R  +  i:4R 

AR2}J  =  4^25. 
152  AR^i^a-i^-ar M 

^Kf2=AS2K-^RP^ 
C    CHECK    FCP-     LL*cr.    CALl  cSATlTN    RA.Nofc 

IF  (4  9"i.L'l   .^L  )i      L  1     TO     lbO 

IF{«ih^.LT.£r«LCJ     Ul    TO     IbO 

K=  1 

P  E  -I  ET  R  AT  I 


Z2  (40  f<iO/  ,X3  (40)  ,Y3  (4U  J  , 


CALCULATES  1NITI'-L 
CALL  rENFTSUl 
RAJL=A4RE3  {  1) 

RA,M^kc.  (2  ) 
CALL     PENPTS(..l 
P3JL=^^ntS  {  1  ) 
i2  J 


:n    POINTS 
pf\2  ,<1   ,Y1  ,Z1  ,Hrl  ,ARn  ,£R?ES  J 


*2  iA2  ,Y2,Z2  ,HT2  ,-iKM2  ,ARREC.) 


^,X2  ,Y2,Z  2  , HT  2, A  K.12,  A  R  PES) 


164 


RB  j^  =  £*R  t; 
CALCULATED    ^LCCE^IVt     °ENETRATICN    PCI  NTS 

CO    140     I=1,1UUJ0 

AR2P  =  ^A2rt  + JA- 

ARM1=AR2R-*R  PI 

0RN2  =  £R2I<-AkP2 
CHECK    FC«     LPPE>v    laH  tP^TICN    RANs»E 

IFIAPMI-GT.A^UPJ     GT     70     195 

IF  (  -l.hN2.u~.  i.-UF)      U     TO     190 

CALL    PE.\PTG(.NL,.N^,Al  ,Y1,Z1  triTl  tAKMl  f ARRcS) 

PAlL=£P^tS  I  1) 

PAiR=ARRtSI2j 

CALL     P-UP  i  »(  41  , 

R8IL -£JRE"  ( 1  I 

R81R-ARRESI2J 

IF(RA0L  .L  J  .  100J 

IF  (R tJL. tC. 

IF( k3Ja.EC. 

IF  (RAJR.t'l. 

IFiRML.EC. 

1F(RFIL  .E  3. 

IF(FEU.t). 

I F  ( R  A  I  * .  E  £ . 
CHECK     F03    JlREr, 

IF  US  L.  M:.  10) 

ISL-  I 

GC    TC    idJ 

I  F  A  F  1 1  L  .  N  l  .  R  H  I  L  ) 

RESL(K)*>~~iL 

PES2  (K  »-  *r 


10J  ) . 

1  ^  L'  J  . 

1  )  00. 

l  o  o  u .  ) 

i  J JO  .  i 

1)0).  i 

T    HIT 

n 


I 

J   t 
T  J 


I  3L=  10 

IS  L=1J 
I SR  =  lU 
UP=  lu 
i"L  =  10 
I SL  =  10 
lSR  =  U 
R  =  1C 
LcFT 
164 


INTERSECTION, EVALUATE     IF     RECUlRfcD 


TC    160 


IF (K.fcO.2) 
K=2 

GO    TC    13  j 

Check  f  jp   i  nz. 


J    TO     155 


"hJlATE     LEFT     I  MT  L  RS  tCT  ION  ,  EVALU  AT  E     IF    REQUIRED 


180  IF(  ISL  .  it.  io)     GO     T  j     JfcJ 

IS  1  =  1 

GO    TO     i  if  i 
160  IFUFAJL-RBJL  )/  (  :  '-,  IL—R-  IL  )  .3T.C.0) 

*A=l*r2r  J-^2  *)/  iPAJL-RA  1L  ) 

*B=(^P2PJ-Ar,2'  )/  (^MJL-RRIL  ) 

B0  =  AR2*-M-i-R  i  1  L 

RES1  (K  )=('<<-  <  \)  /{  AA-AS) 

FES2  (K)  =  J<i*RES  l(h  J*E/: 


GO    TO     ldJ 


101  C2U0 

[NTC2  120 

[NTc2  1-;0 

INT  02  140 

IN  Tu2  ll>0 

N  T  C  2  1 1 0 

NT021  /o 

'  J  T  C  2  1  F  0 

IM  1021  '-.0 

[NT022  00 

NTC2<:  10 

I N  T  j  ^  2  2  J 

NT022jO 

NT  J2t*J 

'.  NT  U22  5J 

IN  Tl^£0 

[NTC^2  70 

;NT022^J 

[  N  T  C  d  2  c.  C 

INT J2i 00 

NT  023  10 

[NT  C*  2  20 

[NTo23:>J 

:NTC2i'.C 

INTC23aO 

!NTC2  3'.>0 

NT^i  70 

INT023OO 

INT  02 390 

[NT02tO0 

NT024  10 

.  -M  T  C  2  4  2  0 

[NTC2420 

NT0244J 

[fJTC245C 

IN  T  C  2 1  6  0 

[N  TC2^»  7  J 

IN  TO 2 4  6  0 

NT02490 

NT02  50C 

NTO^b  10 

NT02520 

NTC2  5  30 

NTC25hO 

NT 025 5 J 

NTC2560 

NT  02b  70 

NT025yO 

NTC2t)S0 

NTC2600 

NTC2C1C 

NT02620 

NT 026 j J 

N  T  C  2 1  4  0 

NTU26  5J 

NT 02 6 60 

NT  0  2c  .'0 

NT  02  6  00 

NT0269J 

N  T  0  2  7  C  0 

NTJ2710 

N  TC2  720 

N  V  0  ^  7  20 

NT0«i  740 

N  TC2750 

'\IT02  7t>0 

NT02  7  7O 

MTO<:  (BO 

NT02790 

TC28C0 


[N 
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FILE 


VOR 


183 


IFtK.EO.^)    ,i 

K=2 
CHbC  *    F  IP    ui  PctT 


;  tj 


B   1 

1^5 


Nib.  10) 


hlT 


It  I 


182 


162 


IF  (  I  jrf. 

GJ  TC  L  f-  J 
I  r  (  R  A  H  .  j  E  .  R 
REr.l  (K  J=^d  iR 
RE  <2  (M  =  As2  < 
]F(K  .fc  J  .* )  tj 
K=2 

GJ    TC    16  J 
CHECK    F1S    i  ITc  HcilIITE 


C  J 
7  J 


"IGmT 
161 


IMT  ERSECT  I  IN,  EVALJATE     IF     REJ'IIRED 


Us)     GJ     TJ     182 


J    TJ     la: 


P.  M:.  1  J) 


T  1 


1GHT 
162 


IF  II 
IS3=  1 

GC    TC    ibj 
IFUP*»JF-PcJ.<)/(F 
«A={AP2:?J-',.K2R  )/( 

;>a=(/>F2Rj-;u2k)/ 1 

pb=4R2<-a-t*<-'  i ; 

FESl  IK)  =  (Ei-eAJ/  ItA-tiB  ) 
RE  S2( K) =Ah*«l SI I n) f dA 
IFU.E3.*)    tP    TJ     155 
«  =  2 

C     RE -INI TI ATE 

163  bHtf-J=  (■*** 

PAjL  =  F.AU 
RAJR  =  MIR 

pe.iL=^d  il 

14C  R6j:<=P-ilR 

KR  I  T  E  (  o  ,  2  i.  J  ) 
210         FQf^ATC     l:CP 

GJ    TC     15b 
19  5  IF{«ESK  J)  .E'J 

IF  (ntSl  (2)  .tC 

IKRES21  ll.cJ 

1F(RES2(^).EQ 

GO     TL     155 
C     £y/ALUATIfjl.     "IF 


INTERSECTION, EVALUATE     IF     RECUIREJ 


M^-SRIP  )  .GT  .O.J  ) 
M^-«  Is  ) 

?H  JR-Rol  R) 


GC    Tl     lo3 


14J    T 


•  9  0v.CC) 
.  JOOJ.O  ) 
.  SClC.  C) 
.90  J  J.  J  ) 


hRCUGHi 

GO 


."-JO    PC  INT S     FOUND1  ) 


\1  L. 


i 


71 
T  C 

ra 

n 


191 
191 
19  1 
19  1 


1  JTER  SfcCri  JOS    WHEN    CLRVtS    ARE    ALMOST    TANGEf^T 


i'-Ut-\  If  A3,  Yj,Z  JtHTi,VGn,APRE  S) 


19  1  YGC  =0.0 

IEPS-=2 

CALL    FENP7 

Pi  =  /  ES  (  *  F  R  t  S  (2  )  -  A  F  R  E*  (1  J  )  /2 

CALL    Pt  19  I  S(  I2»N  1  ,XA  ,V-.  ,ZA  ,rtT2»YGO,  ARRcS) 

R2=*e?  {4.H  ES  U)-ARRES(  ill/i 

RE  52  11)  =(  (  tU**2)-(R2**2  )  +  (  AriP2**2)  )  /  (2.0*AR=>2J 

RE32I2 )  =  RE52(  1 ) 

PtSl  (li  *Su'kT  UeS  U-<l**2)-(  ?ES2  (1  )**2  ))) 

RESK  2)  s-kc  SI  I  I) 
155  PETIJRN 

ENJ 

SUBs  lUTME    PENPTSI  I  X  ,  J  Y  ,  X 

01  "L^ii'J    i  (tjj  ,  /.M  <t^)  ,Z  { 


IR  =  l 

C    SEARCH    FOP    J     JF    LOWER     r    LIME 
OC    101     1  =1, \Y 
IF     (Y(l).Gr.YwJ      >C     TC     102 

101  CCNTIf'Jc 

102  d=l-l 

C    COARSE     jCAI    F.1R    ZERO     °\S5 

C    NEXT    4    LT'ci     ARE    EXECOTEJ    II    FIRST 
>C  =  {  Yb-Y(J)  )  /(Y(  j  +  1  )-Y  (  J)  ) 
ZX(1  )  =  UI  l,J+i)-Z(  i,  J)  )*YC+2I  1 
IF  ilX  U)  .fcC.hT  J     CC    T  C     U7 
13=2 

112  OJ     lu     1=13, 4X 

ZX  (  i  ,  =(Zt  I  ,  J*l  )-<.  ii,  J)  )*YC+Z( 


,Y,Z  ,HT  ,YG  ,XR) 

♦  0,  *0  ),  XI  Aj)  ,XR{  U) 


SEARCH    INLY 
,J) 


Z"4  =  (  Z  I  I  ,  J  «■  i )  -  Z  (  I  -  1 ,  j  )  )  *  (  Y. 
IF  <Zx(I  )  .E'T.nf )     lJ    T  J     lu8 


i.J  ) 
-Y(J))/lY(J+l)-Y(J)i+Ztl-l,J) 
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fILt:    VJP 


FORTRAN       PI 


N 


110 
111 
113 


.0    TQ    114 


GG    TC    HI 


,Y,  Z.X^  ,YG  ,\X,NYi  ZRES> 
ULib-UTllN    SU'F/.CE    TH    EVALUATE 

ULT    ^ETURNtO    IS     ZRES 


IF  ((ZX(  1-1)  -HI)  /(  ZX(I)-h7)  .LT^Uivtf)     GC    T 1    104 

10J       continue 

jC    TC    llo 
C    CALCULATION    UF    PENc  tf  ATI  GNS    IN    EITHER    TRIANGLE 

io<»       x*=x<i-iJMxi  d-a(  i-ii)*ho-tui)/ivijtii-nj) ) 

IFlZN.tC.bT)     «C     TC    105 

1F((  Z  X(  I  )-n  n  /(  /  ,i-hTI.lT.u.U)    OG    7C     1C5 
C,Ll  =  ^tc  U-r-L^+.',^(  ZX  l  I-,.  I-HT  ) 
)S=MI-1  }  +  (  X  -•,-.<(  i-1  )  )*Al>SI  ZX(  I  -1)  -HTJ/OD 
GC    TC    i 10 
105  XS  =  XV+  iX  u  J-A^U^dS  (ZV    -hT  )/  (  APSIZ'-ETl  +  ARSIZXI  I)-hT  1) 

gj   in    lie 

C    LOGICAL     FtNbTR  AT10NS     ACCUMULATION 
1C7  >K(l)=XHJ 

1=  1 

GG    TC    11 J 
lOfc  >R  4 1  R  J  =-  x  (  ij 

GO     11    111 

X^  t  1R  J^XS 

IF( 1R.E0.2) 

IR  =  2 

ie=i+l 

Gn    rn   n< 

116  IE  (1R  .tc.i J 
XR  (21=1000.0 
GO    TC    xln 

117  XR(1J=1JJJ.U 
*-<(  2)  =  Lk.CC,  i 

114  PETUPh 

ENG 

SUBROUTINE  IslTPLlIX 
C  LINEAR  1MERPLL/S*  IDN  CM 
C    CP     VALLE    AT    Cl-ukJINmTcS     XGiY^.THE    » E 

GiVtNSUN    AhJI,VCiu),ZUo,4C) 
C    SEARCn    ECR    I     CE    LcFT    X    LINE 

00     1    1 C  =  1 1  N  X 

IF  Ix(IC).CE.aG)     GO    TO    2 

1  CCNT1NLE 

2  I=1C-1 
C    St*PO     fCR     J    CE     LCWEF    /     LINE 

00    3    JC  =  1  ,  ii 

IF(Y( JC  l.oC.Yal     GO    TO     4 

3  CONTINUE 

4  J=JC-1 
«L=(Y(J  ♦!  J-Y(  J)  J/(X(  1*1  J-X  C  I  )  J 
BL=Y  (J) -^l*X  (1) 
YCR=<AL*>G>  ML 
IF  (YG.EG.YCK)     GO     T)     7 
IFI  Yj.  j  T.  YC-O    ou    TG    5 

C    RESULTING    Z     13     I  4    LJkRER    TRIANGLE 

XJ=XUHJ 

Y3=Y< J) 

23=Z  (1*1,  J  ) 

GO    TC    b 
C    RASULl'INl    Z     IS 

5  X3=X(1J 
Y3  =  Y(  Jfl  1 
Z3  =  Z(  1, J -HI 

C    GENERAL    CALCULAT1CN 

6  Al=/(I) 

Y1  =  Y(J  i 
Zl-Z  ( 1  »  J) 

X2  =  X(  1*  1) 
Y2  =  Y  iJtl  i 
U^lk  1  +  1  .  J*l  J 

AP=(  (Z  1-ZJ)*(  Yl-Y  2)-(Z  1-Z  2  )*(  Y1-Y3J  )/(  (  X1-a3) 
1*(Y1-Y3  J  J 
6P  =  (  <  Z1-Z2)  --«P*  (  X1-X2)  J/IY1-Y2  ) 

CP  =  Z1-4P*A  i-  J,r>*K'   1 

ZRES  =  £P*X^*-EH*r  L  +  CP 


IN    UPPER    TRIANGLE 


LCCC    FCR     eOTE    TRIANGLES 


m-Y2)-(Xl-X2  J 


NT C 3  bio 

NT03520 

:  N  T  C  3  b  3  C 

NTGJb^O 

MTGJbbO 

IN  i  C3ocG 

.NTJ  jS7J 

IN  TO  3b  EC 

INTO  3590 

NTG3oJ0 

N  V  C  3  6  1 0 

NTJ  io2  > 

[NT0J63J 

[N  TC34.4C 

NT  J  Job  J 

,NTC3t'6J 

[NTC3c7u 

,NT0  3o30 

N7Q3690 

NT  03  /UO 

NT0j7  10 

INTCJ720 

INT03730 

.NTC3  74  0 

NT  0  37  tO 

;ntq3  /oo 

4  TC  37  70 
IN  TJ-J73J 

NTG37'^0 

N  T  0  3  8  GO 

[NTOjbl^ 

MTC3S20 

[Nf  03830 

INTC  3840 

[NTC365C 

INT  03  860 

NTG3870 

[NTLifcEO 

[N70  3G90 

NTC3SCG 

NT  0  39  10 

NTCj920 

MTCjS30 

\JTC3940 

i  MTU  39  50 

NIC  3900 

NTJ3970 

N  70  29  30 

[NVC3SS0 

.  M 704 CO J 

1 N  T  C  4  C 1 0 

IN  rG<t020 

NTU4G30 

[NTC4G40 

NT040  5J 

.NTJ406J 

:nt  0^0  70 

NTO'-tOBO 

NTC4CS0 

INTC41C0 

[NTO^l  lu 

N  T  C  k  1  2  C 

NT  04130 

N  7041'*  J 

[N  TC4150 

.M  l"04loJ 

NT04170 

:ntoi  ico 

NY  0^1 vO 
.  N  T  C  4  2  C  J 
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FILL:     VJ?  FORTkA'J       °l  N  P 

»!t|«IV 

GJ    7C    IOC 
C    RESULT  1  nIu-    Z     IS      U     3IVI.MM0    LMC 
7  4*bS*(iii»LtJ«-lJ-4UfJn*(Xu-XUn/(X(I*l)-X(I))+Z(IfJ) 

IJC         CE  fjRN 

£NC 

SLERJIT  iMt_     N^XARY  (Z  t  2  "AX  .  MX  ,  NY  ) 
C    EVALUATION     IF     MxUlH    CR    Jf    CALlB^aTIJN    SLRFACl     ^aTRIX 

03     11  =  1  ,  JX 
OQ    1    J-l,  JY 
i  IF  (HI  ,  J).  oL.Z^mA  i     Z"AX=ZII,J) 

RETURN 

ENC 

Sit* 'JL  TINE     ""UNmSY  (Z  ,ZMNtNX  ,  NY  ) 
EV4LU4ri)M    OF     M  I  vj  IM  J  ;«•    "  P    O    CALIBRATION     SURF  AC  fc     U  TR  I  X 


Q'J     1     I-l,i\X 

go   i  j=  i,  i> 

IF l/ll , J) .Lc Z*1M 
RE  FlMN 
ENC 
END    CF    VL'R 


Z-^i-J=  Z  (  i  *  J  1 


INTO 
INTO 
INTC 
INTO 
INTC 
I  NTC 
INTO 
INTC 
I N  T  0 
INTC 

rue 

INTO 
INTO 
INTC 
I  MY  J 

INTC 
1 N  i  0 
INTO 
INTC 
IN  T  (J 
INTC 
INTC 
1NTJ 


4210 
42  20 

42  iC 
4240 

4  <:  1 0 
42  70 

42  EC 

4  2  '*  J 

tOOO 
4310 

♦  3  20 

43  30 
-.3  40 
4  3  jO 
42C0 
4  3  7  J 
43  JJ 

43  SO 

♦  400 

44  10 
4420 

♦  4  JO 
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APPENDIX  2 


LISTING  OF  SWVDR 
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t  I  / Ct  /  <  S        it.  4  5.  19 


F1U:    S^VCR 


F CRT  SAN       Pi 


N 


C     JTART    3F     S*vOR  ••»*•• 

C  1  '^  fc  N  S  1  1  J    Alt  (  lu),jhll  10)  ,  >A(  4C)  ,  YA(  4C)  t 
iZA(<»0l'*0),<dt4JiffJ(tJJt/i(t0,4jj,Xf«X(^J),rAX('*J)i 
2ZAX  4o, 4C)  ,  M  a(  -h.)  .YlM^tUl   ,*SX(40  ,-tOJ 

119  FCF"^TCZli>i 

131  PurxV^Tj  f  7.^  j 

13 C         f-THlAT(F£.  1) 

120  FChV£~  [  fd.'j  J 

132  KH"MIF  1-).  si 

135  FCF'-1AT(  Fli.o) 

136  FGH."*T<FW.S) 

C    KEAO    P°  -ME    C^L1^?ATI-j4    DATA 
Rt*C(2,ii:»)    .  I  At  4  Y 
ON     121     i=i  ,  IX 
Rt'AD(2,  1JJ)     XA(  I  J 

121  Ci-NTI.VJE 

0J  1**  I  =1 ,N  r 
RfcADU,  lij J     rA{  I  ) 

122  CCMlKJt 

DJ     1*2    1=1  ,WX 
CO    123    J=l,  ii 
Ffc*D  12,120)     l*(  I  ,jj 
122         CGN'flNLc 

CO  221  1=1  .nx 
REAJ( 2,  120)     Xd(  1  ) 

221  CONTINUE 

ZQ  222  1  =  1  ,NY 
READ(2,  13CJ  YB(1  J 

222  C'JNT  IfU  E 

DC  22'i  1=1  ,  NX 
03  lt±  J=l,  JY 
FEAl  (*,  UOi     ZE(  I, J) 

223  CUMlME 

C    TRANiFORV    C^LInR'WlJH     MATRICES 
CO    2  2<*     1=1, NY 

224  XAMI)=VMi) 
DO    *:*5    1  =  1,  NX 

225  YAX(I)=X*U) 
OJ    226    1=1, NY 
DO    <r*t    J=],  Na 

22  6  ZAX1I , J) =  2MJ,1 ) 

J  J    **7    I  =1  ,.\r 
227  X6X(1)=  Yfc  (  1  ) 

CQ    22b    i  =  l  ,NX 
22  E  Yd  >(  1  )=>P(  I  ) 

CO    2  29     1  =  1,  J  Y 

DO    229    J  =  l  ,mX 
22  5  Z9X(  I,J)=Z?(  J,  I  J 

C     FRCBE     StTTlM£5 

YPA=C.O 

YP8=55. 

XPA=O.J 

XPC=2C. 

XPC=-20. 
C     FKCOPAf*    CCNSTA.NTS 

Y"LO=-JC. 

YRIJP  =  3j  . 

Yk1N=-30. 

XRLO=-3  0. 

XRlP  =  _iO. 

><  IN =-3  0. 

ICHPS=  1 

NCCCFS=2 

"M0S1M=2 

IIT=1 

ISCAN=1 

REl>PS=0.3 
C     EXPERIMENT     ^  I-1JLAT  ION 
153         C'INllNLc 
401  FJP^ATf F15.H) 


[NTOuOlO 

:ntojo^,) 

i  M  T  C  C  C  3  C 

INT  J  J  J/.O 

NTOOO^>0 

N  T  C  0  0  ii  J 

NTJ0J70 

.NTCOChO 

NTO  JJ^JJ 

NT  J  j  I  00 

[NTOOl  10 

[N TOO  120 

NTJoL  jj 

INTO  J  140 

NTJ013U 

M COl  tO 

NT  wC  1  70 

;urooidO 
intcciso 

IN  f^02o0 
NT00210 
[N  7<.C**0 
INTO  02  JO 
INT0O240 
INTOJ250 
NT  002 60 
;  N  T  0  0  2  /  0 
[NT00260 
:NTJ0290 
[NTCC3CC 
[NTO0310 
.NTJ032U 
[NTCOJjO 
;NTJu34j 
;NTCC35C 

;ntoojcO 

NT  00  3  70 

NT002&0 

IMTCC39J 

.NTCC41C 
[NTu0420 
;NT00430 
[-1TCC440 

1NT0G450 

;ntoo4oo 

NT0O47O 

Ni'0u4  fcO 
;MTCu4  90 
[NTC05C0 

;ntoj51o 

N7C^520 
NT0JS2O 

NTUJ5 V) 
>J  TCG5  50 

[N  "Oo5oO 
NTOJS  fJ 

'.  NT  0^5  80 
NT  JO 5 90 

;ntco6cc 

[NT0C610 
NT  0 jo  20 
N  T  C  C 1 3  C 

INTO  )6  40 
N70uo50 

:^itcl.c60 
;ntooo/o 

NT006  30 
[NTG0o90 

ntoo/oo 
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FILE  :    5WVDK 


FCRT^mN       PI 


READ(5,<tLl)     PT 

WRlTE(ciJl)    ST 

IF  L  ^  CS  i  f.tw.2  J     oC    TO    500 

PEAOISi,  «tl)     PS 

WtUTE(6,  tJi.)     PS 

Rt<0  {5,  ■■♦CI  t      \LI"A 

write(c,ald   uf'; 

PEAC(5,40i  I    kHI 
*RI  lfc(  c  ,401 J     Prii  i 
CALL     !*JTFLT(X«f  Yd 
ALFC=;LF*-aPC 
CALL     hNlt-LT  1X4,  Y^ 
/Lfl>£L  M-XPC 
CALL    I  MfLT  Ia<s,Y 
PF18=PHl  I- YP3 

Call   intclt  uh,y 

poyn=pt-ps 

ppa=cpa*pjyn*ps 

FFe  =  CFj*F^YfwPS 
PPC  =  CPC*PJY  >l*PS 
FFC=CPJ*PCY   U"S 
C     »<EA0    MEASUrcfE-NTS     CA 
IFlNGSiM.t  J.1J      , 

500  PE/i:(5,iJji    3:>m 
*R1TEI»,13M     PPA 
RE/0I5,U.>)    °£M 
WH1T£(J  ,1ju  J     PP3 
READ(5,lia)     PPC 
WR1TE( j,  Jjoi    ?PC 
SEAd(5,l35)     PPu 
WRITE( c,  1  Jc)     PPu 

501  i-tAC  (5,  Up  J    3SI'I 
WRiTti6,l36)     PSiN 
IF(  ISC4N.EQ.  I)    j'l    TO 

180  IF  (iSCAf^.cT  .iOOJ  i    „0 

II  1=1 

ICCPS=1 

1SCM=1SCAN*1 

PSI.\=PaIN+5.  C 

WRIT  EU,  ljfc)    P*il.J 
1  El         CCMlMt 

PS=PSI < 

FSC*PSIN 

PST=PS 

c  calculates  pressure   coefficients 
300       cotik'je 

PO*N=Pf-?S 

CPA=  (PPZ-PS)/PDYf 
CFE=  t  FP  E-F£  J /PCYr 
CPC=(P?C-Pi)  /"JY*. 
CPL=  (PPC- J5  i/^CY'j 
C    CORRECT?    PS    AS£L^FTi<_ls    TC    ENSURE 
IF  (NJCDPS.E  0.2)     00    TO    ->0  1 
2ANpX=-1 JJuJ  . 

CALL    MAXiKY(ZA,Z  AI^X  ,'\IX,NY) 
IF  ICPA.C1"  .ZM  \X/     :5SC=i.J  i.*PS 
IF  (CFC.  ul.  <.*-^x  )     PSC*1.0  1*rS 
lFICPO.OT.Z  '<  *IA  A)     P.SC  =  i.Jl*PS 
2B^X=-iJ  JUJ  . 

CALL    f^AXAnY  (  Zd.ZcMX  ,f\X,NY) 
IF(CP6.GT.Z'3'U,0     PiC  =  1.01*PS 
ZAH  Mi  JJuO. 

CALL     MlMArY  IZ4,Z  AMi,  ,NX,NY) 
IF  (CPA.LT  .It  1  IN  )    »SC-=U.99*PS 
IF  (CFC.LT.ZA'MNJ     PSC=0.S9*PS 


•ft l«i < 


,ZA,ALFA  ,PHII  ,NX,NY,CPA) 

,<1A,ALFC,PHII,!\X,NY,CPC) 

/  Ufi  i*»LFC,  C|-  II, NX,  NY,  CPC) 

F,ZB,ALFA,?H18,NX,  NY,Cor3) 


TO    501 


iai 

TO    ^05 


CF    ARE     IN    C/LiERATilN    naviGE 


J  J 


CC/9*PS 


IF(CPD.L"I.ZA1I 

2B*lN=i JJOO  . 

CALL     MI.MA^Y(  Z-*  ,ZoMN  ,IV<  ,N/ ) 

1FICPP  .LT.Z'VWn     P3C=0.99*PS 

IF  (PSC.  EI..P3T  I     tC    :C    30i 


NTC 

NTO 
NTC 
NT  J 
NT  J 
NTC 
NTO 
MTC 
N  )'  J 
NTC 
N  TC 
NTC 
M7J 
i  TC 
NT  J 
.N  10 
N  I  C 
NTO 
NTC 
NT  J 
NTC 
NTC 
MTU 
NTO 
NTo 
NT  0 
NTC 
NTC 
NTO 
NTC 
NTO 
NTJ 
NTC 
NTO 
NTC 
NTC 
NTu 
NTC 
NTC 
MTU 
NTC 
NTJ 
NTO 
NTC 
NTJ 
NTC 
vlTO 
NTO 
N  TC 
NTO 
NTO 
NTC 
NTJ 
NTC 
NT  J 
MTC 
NTC 
MTU 
N  T  C 
NTC 
NTO 
NTC 
NTJ 
NTJ 
NTC 
NTJ 
NTO 
NTJ 
NTO 
INTC 


C710 
J/ 20 
l  ?3C 
J740 
J?5J 
J7b  J 
0/7J 
0780 
J  l^J 
JrtO  J 

oa  n 

08  2  J 

C  fc  •'•  c 

Jd5J 

000  J 

oe  70 

J3oJ 

casj 
0900 

JS  LJ 
o92J 
J93J 

09  40 
j  V  50 
J960 
CS7J 
CSdJ 
J9QJ 
1CC0 

1J1J 

1  J2J 
1030 
I J  40 
1C5C 

I  Joo 

10  70 
1C8J 
lJ^J 
il  JJ 

11  10 

112J 
1UJ 
114J 

1 15  J 
116C 

II  13 
11BJ 
11SJ 
120) 
1210 
1220 
1230 
1240 
125J 
L/.'oJ 
127J 
12HJ 
12S0 
1  JCU 

iJU 

1.320 
U  J  J 
1j4J 
12  50 
1j60 
137J 
UdJ 
1j9J 
14C0 
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FILL:     SWVLP 


FDRTRAN       °  1 


N 


«'u«i»' 


CF    THE    FO'IR    VALUES 


ICOPS  =  ICCPS*l 

IF  (  ICHPi.GT.  10)     oO    TO    180 

FS=FSC 

PST*PSC 

GO    TC    JJJ 
301         CJM1NU 
C    CALCULATE i     it.  "I    Ai_F     e    PHI     ANtiLcS    F  C-R    •  I*     S    Ml' 

L*LL     INMCS  (  1?  A,  yo  -I,  Y"L  i,  Y  RJP,  YRIM  , 
1CPA,CP^,Al(-  ,  PHI  ,XA  ,Y  A  ,iAt 
2XB,YR,ZJ.  ,A^XtYAA,Z4A,AdX,YHX,ZeX,\X,NY) 

/LFl=ALF  Uj 

AlF2=AlF  (2) 

Fhii  =  PH<  1  J 

FHl2  =  F<ll  (2) 
C    CALCULA1FS     1*3    ALF     €    P,  1 1     A^JuLtS    FOR    M'     R    'ill 

CALL     1NTSCS  UPA.XPCtXRLOrXRUPfXRIN, 
1CPA  ,CPC  ,FHI  ,-ii.F  ,  >n<  ,  Y*X  ,  *.AX  , 
2XAX,  Y.'.A.Z'U.aA,  »AtZA  ,AA,  Y-;  ,ZAtNlf»  >IX) 

ALF3  =  AL»-  I  J.  J 

ALF4=plF 12) 

FHiJ=PHi  (1 J 

FHH=Phl  (2  ) 
C    SELECT^    THE    P<)PfcR    ALr    fi    PHI    ANGLES    CUT 

FALFM  I  ALFi  ♦:LFtHLFj*iLF-t  J- A  UXK  ALF1,  ALF2,  *LF3, 
1AlF4)-av1.iIULF1  ,/:lF2,ALF3  ,ALF<+  ))/2 

FP  Fl=(  (  PF  I  !+►>  rt  I  ^+?-ll  jf!>HI  t)-V'AAl(  Pi1lltPrU2,?HI  i  ,PhI4) 
l-AflMlPm^Fli.Klh.PhhJl/Z 
C    CALCULATti.     JE*    °i,    F«,M    D'\l\     CF    PUjITI]^    0    AN'O    CORRtCTS     P~ 

/SLFCN=r  ALF-X-»0 

CALL     l-JIf-uTtXUfYM.Z^.ALFDi^.FFHI.NXtNY.CPCJ 

AL  FCN«FAlF-XPC 

CALL     1.-ITB1.T  (X  *,f  A,ZA,ALFC"I,FPHI,  >IX,    IY,CPC) 

PSN  =  (P?C*C  PJ-PPJ*CFC) /(CFu-CPCJ 
907         WR  1TE(  a,«,l  1)1  IT,  FALF,F?HI,  P  SN 

911         FUHNflT  I  •     ITfcfuTlL'N    VALUER  •  ,  la,  10X,2F6.2,  FlO.2) 
C    CUUVEkGENCc     fEoTS 

tPSPS=^pHS(  (ajvj-pj  J/PII 

EP  SPSG^.  00000000  CO')'}  <(PT**2.0) -.0)00  O0'33*PT 

IF(  UPoi.t0.2)    cPS»S3».  lCjCuOOCG09*(  FT**2.0)-.  033  0 003* PT 

1FUEPS2.L-C.2)     tP3:,oG=.JUJJOJOJuJ9*(DT**2.0)-.CuOOCC6*t,T 

IF«EPSPS..iT.EPSPSJ)     GC    TC    310 

GO    TO    jjj 

CCM1MJE 

m=n  i+i 

II-  (  I  IT  .IT  .uj    G1    T")     180 

FS=FS+«LLXPS*  4PSN-PS  ) 

PSC=Pb 

FST=PS 

GO    TC    30C 

WR  ITEJo,  Hi.) 

FL->f*iT  {  «0«  ,  ' 

WR1TE(  o  ,  W2) 

FC^ATf     VtLJCiTY 

kFUTtl  6  ,1/J  )     PSN 

FORMAT!  '     nA  f  I  C    PkcSSUkE 

l*KiTE<of17t)     FT 

FJRI-IAK  •     7  JT^L    PsESSlRE 

AM=SCFT  (  (  (PT/PGr 

fcFUTbto  ,1  /?)     AM 

FORMAT!  •     MACrl    i»uN9ER 

CO    T"i    ivj 

WRI  lb (o  ,307) 
-FTRy.  *T{'     jCAWNEU     1J00    TIMES 

fcU"    TC    iOj 

ENO 

SU^RCUT  I  mF     I  JTSC^  (  *kP  1,  ARP  2,ARL0,  \RUP,A^  IN  , 
1HT1  ,hT^  ,  hlS  i,  rtji,,(i  ,Y  i,  Zl  , 
2X2  tVifZ2|Aj,Y3iZ2  fX4f  Y4iZ'»  »'M  tN2  J 

CI^tNS  10  ^     -<L5i(lj),-  E3  2(  10  ),  4 'RE  SI  10)  ,  XI  (40)  ,Yl(iC) 
IZ1  (4U|<tC)  ,  X2  (40)  ,Y2  (40)  ,  Z2  C»0  ,tJ),  X  JUJ  J  ,  Y  3  (  <t0  )  , 


310 


30  3 
171 

172 

173 

174 

J7i 

3C5 

30/ 


FALF 

VlLLCiTY 
FPi-i  i 

P  ITCH 


YAW 
15  ' 

IS* 
IS' 


IS*     ,Fl3.2) 
,F  15.2) 

,F15.3) 
,F1^.3  ) 


)  **0.2t)57l4-  l.C)*5.0) 
IS'     ,F15.5) 

) 
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FILE:    SWVCP  FDRTRAN      Pi  MPS 

21  3(  «,iJ)  ,X4(  <«0)  ,  V+I4CJ  ,Z*(40.»AO)  INTCcllU 

RE31  (  li  =  jjjj  .  lNTJ^l^O 

Rt  <l  (?) =9000.  •  I M  T  C  2  1  3  0 

RfcSi(  1)  =iC0L.  INT  02  1<»U 

PES2  (<f  J=9j00  .  INT0215O 

iSl  =  l  HTC21t0 

IS^=  I  INTU2L7J 

0AP=1.  1NTC21E0 

ARiR  =  A-<IN  INT  0*1  90 

uC    TC    152  INTC2200 

15C  dk?R-=AR2R«-i;AR  HTC2210 

AR2RJ=\k27  IiN'i'0^20 

152  *K'n  =  Ari^S-^5Pl  I-NTO^jO 

ARM2-=A«2R-4R.J2  INTl-2240 

C    CHECK     FOR    L'"J*e.<    rALlBRiTiJN    RANGE  INTC2250 

IM4i^.LT.i)hL:3)     LC    TC    150  INTC2260 

IF(Ar(f2.  L  I.  AriLC)     OC     IC    150  INTC227G 

K=  1  !NTU22dO 

C    CALCULATES     INiTiJL    FtNcTSJTiCN    PCIMS  INTC22SC 

CALL    PLNPTSJ   Jl,u2  ,Xi  ,Y1  ,Z  1  ,HTl  ,AR'-<1  ,ARRES)  INTG2300 

PAJL=ARPES  (  U  INTG23  1J 

RAjR=A*KtS(2)  INTC232G 

CALL    PE.NPTS  (J  i,N*  ,X2,Y2tZ2tHT2tARi-12  t  ARRES)  INTU2330 

KBwL  =  AfhfcS  tl  i  1NTC234C 

RO  jR  =  AkRES(2)  INTG2350 

C    CALCULATES    SUCCtiiVE    °EUETRATION    PCINTS  INTC23oO 

LCJ    14C    1=1,10000  INTC2J7G 

AP^P.  =  AR25  +  J^K  INT02j8J 

^RM1=^k2R-ARP1  INT02J90 

ARV2  =  AK2R-/-rvP2  INTC2<tG0 

C    CHfcCK     FOR     JP^F.R    C  A  L  I  B  R  A  T  I  "J  i  M    RANGE  INT024  10 

IF  (^Pf'i.tT.A-'UP  i     O    TC    19:>  I.NTG2420 

IF(AKf»2.:iT.ARUP)     »L     TG     19b'  INT C 2 4 3  0 

CALL    PL  IPT^l  a»NiiXl,'i  1,2  1,hT1,ARMI  ,ARRES)  INTG244J 

RA1L  =  ^Pc3(1)  -      INTC245C 

RAlP=A.^ESl2)  INT02460 

CALL     PfcNPTS  I'll.  J^iX2fY2,Z  2  ,HT2  ,A  RM2  t  A  RRE  S)  INTC2-W0 

FdlL^ApRES  (1  J  INTC2^fcO 

R8IR=ARRcil 2)  INTu249J 

IF  (P/JL  .E  J. 1000  .  i     ISL=iJ  INT02500 

IF  (REJL.Ei.,1  )00.  )     ISL*1J  INTG2510 

IF(RPJR  .fcQ.  1000.  J     IS?=IJ  INT02520 

IF  (R/SJ-5.EC.  iJJJ  .  J     IS  P=1j  INTC2  5  20 

IF(P.AIL.EC.  10CC.  )     iSL=lG  INTG2540 

IFIRtIL  .fc'Q.  Ijjj.  i     I  isL=  ICI  INTC2553 

IF <RElR.fct.lO  )J.  )     15^=10  IMTC256C 

IF  (RAH.EO.  luuG.  i     IaR=lO  lNT02b70 

C    CHECK    FCP    3IHECT    hIT     11    uEF7     I NTE n SEC T I  ON, E VAL uA TE    IF    REQUIRED  INTC250J 

IF  (i  SL.Nt.l  J)    GJ     TG     I6<t  INTC2t;Sj 

ISL=1  INTG^ouO 

GC    TC    183  INTC261C 

164  IF  (RAI  l.NE.R  U  L)     GO     TC    180  INT02620 

RESKK  J=R-l  »L  INTG2630 

FtS2  <M  =/!r<^h(  IN  Tc2t4G 

lF{K.Ew.2)    GD    Tj     153  INTJ2650 

K=2  I  NT  026  60 

GC     TC    133  IN7C2  6  7G 

C    CHECK.     FOR     iM  IER  MEDIATE     LfcrT    I  NTE  RSECT1   IN  ,t  VA  LU  AT  £     IF     REClIRtU  lNTO^et.0 

180  IF  <1$L  .ME.iO)     G.)    TD     ltJ  '  INT^oSO 

ISL  =  l  INTC27uO 

GO    TO     18  3  INT027  1J 

160  IF  (  (-/ JL-PBJL  )/  (  f-  fiR-RcIL  )  .  GT  .C.JJ     GO    TO     ltt3  I.NTC272U 

AA=(AR2RJ-~R2N/<RAJL-r.AIL)  INTC27  30 

flA=ARtR-AA«- AIL  INT02/^0 

A8=t;F2  ^J-;Mkj/   (CcJL--HIL  J  INTC2750 

BB  =  AR2^-a  3*k->IL  TJT02760 

P£Sl(Ki=(tlu-H4J/(Av-'.rj)  I M  TO.'  7  70 

RE  S2  IK)  =Ad*  PcSl  (I-  jt-Eit  IMTC27EC 

IFIK.EJ.2)    GJ    TJ     155  IN T 02 79 j 

K=2  /  INTC2  8  00 
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FILE:     SWVJR 


FJRTRA^ 


C    CHECr.    FOR    01  REC  7    HI  T 
183  1F<  ISt-  /It.  1J  1    GO 

I  Sh  =  i 

CO    T')    Ui 

161  IF  (  RA1 ?.Nl,=  :3 IR) 
RtSl(tO=tOl* 

RES2(K)  =  \*t* 

if  ih.bw.z  j    g:  to 

K  =  2 

GO    TO    loJ 
C    CHECK    FCF-     lMcH"cCl/iT 

162  IF(  1  cF.NE.  1C)     GG 
IS«=  1 
ou    TC    163 

162  IF  KRAJK-RiJ?)  /(  P 
AA=(  AP«:FJ-AR2R)/  ( 
EA=AR2R-AA*h-4lfc 

AB=(  AR2RJ-A.T**  )  /{RHJP-RR  H  ) 
Ee  =  /SR2F.-£8*RRIR 
Rt  5  I  (  rJ  -  (  ^  i  -  r  A )  /  ( 
RES2  (K  )=a^*RcSl(N 
IF  (K.tG.2i     «C    TO 
K=2 
C    RE-INITIATE 

163  AK2hJ  =  Mr2H 
R4JL=RAIL 
FAJR=PAIR 
PaJL=R61L 
R3JR=P^  l.R 

140  C  INT  ir.'iJ  £ 

t»Kl  7t(6  ,210) 
210  FO^-iTI '    LJO°    IhO 

GC    TC    la  5 
195  IF(KCSU  1)  .Ew.iCC 

IF  (K  ES  1(  <>  .  EO  .Vwi 

IF1PES21  l)  .l'0.9JG 

lMRc$2(  I  )  .fctf .*0u 

GO    TC    15  3 
C    EvALLATI^K    GF     IMERSE 
191         YG>U.J 

CALL     FcNFTS  IN2tNl 

R1=ABM  AR-.c3(2)  -£ 

CALL     FENPTS  (M^tM  i 

R2^*ES{  Jt<RtSl2)-- 

RES2(1)=<  (-  i***)- 

RES2  (2)  =  FES2  U  » 

RE  <1(  II  =  S>JkT(A3S  I 

RESl(2)  =  -RcSl<l) 
155    PtTUfrN 

END 

SDORCUTI  JE    PENPTj 

oi  Ntrs  icr\  y  c+oj ,  * 

IR=1 
C    SEARCH    Ftf.    j    .if-    LOWER 

DO    I  CI     1=1, NY 

IF    hUJ  .Gt  .YG)     o 
101         CONTINUE 
1C2  J-I-l 

C    COAPSc     SCa  I     FIR     ZERO 
C     r>EXT    t     Li  N't1?     f!kfc     tX  t 

YC=(  K—  Y(  J)  ) /I Y( J 

2X  (D-  (Z  (  I*  J+  ii-Z 

IF  UM1I  .  tC.hT)     o 

13=2 
112  CC    ICs     1=  IE,   iX 

zx(  n-(  zi  i  ,  jmi-z 

Z<=IZIl,Ju|-/ll- 
IF  (ZX(  I  )  .tG.hT)     o 
IF  (  (ZaU-I)-HT)  /( 
10  3         CCNT1MJE 


!H     RIGHT     INTfiifJi-^TlCNttVALJAf  i:     IF     HEC'JIkEG 
T  )     ifcl 


C-0    TO     162 

1J3 


E    RIGHT     INTERSECT  1~)N, EVALUATE     IF    REUUlRcD 
TJ    lo2 


UP-R3IR).GT.0.0)     GC   TO    163 
RAJ*-}  A  U  ) 


1  +  3  *V 
155 


FHRJ'JoHtMT    POINTS    FOUND*) 

CGI  oO     TC  191 

C.J)  G.I     TO  19  1 

J.J  )  u]   TC  19  1 

J.OJ  ^J     i'O  1SJ  1 

CTICNS    WHEN    CURVES     /HE    /iL^IST    TANGENT 

,X3,Y3,  Z3  ihTlfYCH, ARRES  ) 

??. E  ;(l))/2 

,X4,Y4,Z'+  ,HT^,YGD,ARRES) 

RscS(l I )/2 

{-<*:**2)  +  l  ARP2**2))  /(2.0*AR?2) 

lRi**2)-( RES2 (1 )**2 ))) 


(  sIX, NY, X,  Y,  Z,NT,YG  ,XR) 
X(4J)fZUOf<*JJ,X(40),XR(lj) 

Y     LINE 

J    TD     102 


p\SG 

JTtC     1M    FIRST    StARCF    OMLY 

t-L)-Y(  JJ  I 

( i,  j  i  ;*vr.  +zi  oJ  J 

:    T1    10  7 


UtJ))*YC  +  ZUtJJ 

1,  J  J)*{  YG-Y(J))  /(  YIJ+l)  -Y(  J)  )  +ZU-1  ,J) 

f    TC     1J8 

ZX(  I  )-hT)  .Lr.O.  J)     GC    TO    1J  ♦ 
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FILE:    SWVJR 


FORTRAN       PI 


104 


GO    TO    lit 
CALCULATE    JF 


PFNET3  Af  IO.NS     Im 


E I IHCR 


TRIANGLE 


XP=X(i-l)MXU)-XU-ii)*JYG-Y(JJJ/mJ*l}-Y(J))       • 
IF(ZM.EJ.HD    G  J    T)     L  C  5 

1FU2XII  )-nT)/(Z'i-FT  l.LT  .j.j)    GO    TO     1C5 
GG=At£  Ihl-ZV)  «acS  (ZX  i  I  — L  J-HT  J 
XS=>(  1-1  )♦(  X'1-X(  i-U  )*A  3 o(  ZX(  I-1J-HT)  /JO 
CO    TC    UJ 

>S=>^(X(l)-X^)*AdSlZN    -MT  ]/(AtiS(ZM-hT)  +  /eS(ZX  m-Fl*  1  J 
UC     TC     iiJ 
C    LClWCaL     FLlvLT^TiCNS    fi  CCL  V'JLAT  1CN 


1C5 


lu? 


M  I) 


lot 

11C 
111 
113 


116 

117 
114 


GO    TO     11< 


GO   TC    117 


>*(1)' 

1=1 

iiiJ    70    113 

*r  (  ir  )  =  >(  n 

GO    TG    ill 
>R(  I  R)=XS 
IF(  15  .EQ.2) 
IK  =2 
IB=I+1 
GC    TG    ii^ 
IF  I1R.  ten 
XR(2)  =  1CCC.G 
GO    TO     114 
>K(1 )=10J0.3 

AR  ( I  \  -  I  c  c  c .  c 

PETUFN 
ENC 

5URRCLI  INF     I  JTPLT(X  ,  Y  .Z  t  XG  ,  YG  ,  fs  X,  N  Y  ,  *  Re  S) 
LINEAR     INTcRFCLATIC*    CN    CALIF-  *T  I  IN    SURFACE    TO    E^LO-xTE 
:P     VALLE    /T    CO-'JKOl.NAfti    /'j.Vj.     TnE    RESULT    PETUR^tG     IS     ZRES 

OlVfcNS  ICN    X  (<«u)  ,  Y  Uu  )  ,Z(  tu  ,40) 
SEARCH    FCh     1    CF    LtFT     X    LINc 

JO     1    IC=i,NX 

IF  (XI  IC  )  .Ll.XGJ     GO    TO    2 

CCMIME 


1=  IC-  1 

SEARCH    FCP 


oo  j   j:=i,ky 


J     "IF    LOWEP    Y     LINE 


YGJ     GO    TC    <\ 


IN  UPkER  TRIANGLE 


IF(Y  (JC  J  .oE 

CCM"  iMt 

J=JC-1 

AL=(V(J  +  l)-r'(JJJ/(X(  in)-Al  I  )  ) 

EL=Y  t  Ji  -AL*X  (  l  j 

KR=(AL*/J  }*HL 

IF  (YC.EQ.f  lR  J     G3     T)     7 

IF ( Vo.^l.y .R)     GJ    T]     5 
RESULTING    Z     IS     I  J    LQURfcR    TRIANGLE 

XJ^X  IIU  ) 

Y3  =  Y( J) 

13= I {  I  ♦  1 , J  ) 

GC    TC    6 
RASLLTINo    l     IS 

XJ=X  (1  J 

Y3  =  Y ( J+l) 

Z3=Z(  1,3  +  1) 
GENERAL    ULCUIAT  ICN:GClD 

>1=X(I  ) 

Y1=Y(J) 

2l  =  Z  11  t  J  J 

X2=x(  IHI 

Y2=Y  (J*  1) 

iz  =  m  +  i  ,.j  +  i  j 

-AP=(  (Z  1-Z  JJ*(  Yl-Y2)-(Zl-Z2  )*(  YL-Yj  J  )  /  {  {  <  1-X3  )  *  ( Y  I  -Y2  J-UI-X2  ) 
l*(Yl-Y3i ) 
EP  =  UZl-Z2  )-AP*(Xl-X2)  )/(Yl-Y2J 

CP=Z1-AP*X  l--iP*Y  1 
ZK  ES=AH*Xi;+rJ.>;*Y^+C  > 
GO     TO     100 
RESULTING    Z     IS    0  1    JIs/UING    LINE 


FOR     eOT  F    TR  1  ANGLES 


[NT  C  3  tlG 

.NTO  J520 

N  T  C  J  5  J  C 

[NTG3  5'tO 

NTC3L>50 

NTCJ5cO 

fjTOJS  TO 

;NTCJ5cO 

NTC35?0 

NT0J6JJ 

[N  TC3t  10 

INTO  36  20 

r;  TO  lb  JO 

;nt  dio^o 

NT  J  in  50 
.MT0366J 
INTliJc  10 
'.  NTGJbbO 

NTC36S0 
IN  f u  2  /JO 

:ntoj7u 

NTC2720 

INT 03  7;>0 

N  T  0  3  /  40 

.  N  T  C  J  7  5  0 

IMTOJ  7t.J 

.NT0J770 

[NTC3760 

NT03  7  90 

N  T  C  3  E  C  C 

lT03d  10 

N7C->d20 

NTCiEiO 

iUTJ3  8VJ 

[NTOJiibJ 

N  T  C  3  1 6  0 

NT0io7J 

N  TCj«3 J 

NTCi'iSO 

NTO3V00 

i  1 1  T  C  t>  S  I  0 

INi  02  J  20 

NTCj^iJ 

[N  TL3940 

[NT  039  53 

NT0j960 

..\TUj9  70 

NT0J9  60 

NTCJS90 

NTCioCO 

:NTj40  1J 

N  T  C  h  C  2  C 

N  T  0  4  0  2  J 

NTCJ4o40 

INTClCSJ 

,  NTO-VOo  J 

NTC4C7C 

NT04O90 

NT  J40-V0 

NTC-ilCd 

INTC411J 

NTO^tl^O 

NTC413C 

;nto414j 

NTCtlSJ 

r  c  h  l  l  o 
;nto417o 

NTC41EC 

iNT0^l9J 
NT  34^00 
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FILE:     SWVOR 


FORTRAN       Rl 


7  ZREs  =  <zu*i,j*n-m,jn*  Uo-XTn  )/(x(i«ij-x(iji«-^<i*ji 

100  7fcTUR* 

t-r\o 

!>U*R  30TI  NE     MAAAn,  Y  (Z  t  Z  -lAX  ,  M  X  ,  NY  ) 
C    fcVALUATICN    °F      'iXMUM     :P    JF    CiLI^KATUN     SURFACE     MATRIX 

Oi^LNSICN     il^.'in) 

DO     1     1=1,  J  X 

CD  I  J* If NY 
i      IF  (  Z  <  I  ,  J)  .  jt.iMAXJ  ZNAX-Z(IfJ) 

RETURN 

ENL 


iL8P.ru  J  I  I^C     Vii\.AhY(Z,2MIN,MX,\Y) 
EVALUATION    IF    'MNHijm    rp     }F    CALIhPATION 

Ji  ^t ^Si  cn   *  hJ  ,40) 

DO     J     I-  I ,  N  A 

CO    i    J=it  JY 

1F(Z(  It  JJ.Lc.^yiNj     ZMN  =  Z(J,Ji 

RETURN 

£I\C 

END    OF     Sk»\*OR 


SURFACE    MATRIX 


INTC 
INTO 
INTC 
HTO 
INTv. 
IN  TC 
INTJ 
INTO 
I N  i  J 
I  NTO 
INTC 
I  NYC 
INTO 
INTC 
INTO 
IN  TO 
INTl 
INTJ 
INTC 
INTO 
INTO 


12  10 
•♦^O 
<i<L  JC 

t^ij 
4ct>0 
42t0 
<*2  7  J 
♦  2  30 

•'♦3  00 
lilO 
•ij^O 
O  JO 
4340 
4ji>J 
-t  360 
s2  70 
4  380 
i3c.C 
-.400 
4410 
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APPENDIX  3 


LISTING  OF  WVDR 


L7/Lt/iS       li.^ii.il 


FI Lfc  :     ^Vl« 


FCRThM       Fl 


'J 


C     STAP1    JF     *\, 

C  i  v£  ^S  1 

l/A  1  +  0  f  ♦ 

2ZA  A(  HiiU 
119  K'RMT  <2  J  >  < 

13  1  nr>*IMIf/.ll 

UO  F'JiJv^TiFJ.s) 

132  F.J<Mil(r  1C.  3) 

13  5  FJR,'*T(Fl3.5) 

13  c  FC*MT<»-l/.bJ 

c  re 40  ° :  i i-t   '^liv.hi 


-*  I'll*!' 

4     iLrllJlf'HK  1 J »  f  XA  (  4  U  t  Y<  (AC)  , 

J  ,  <  d  14  J  i  ,  Y  J  ii  J  ),  ZtM-tO,  10  )  iX  Aa(  tO  } ,  r 


MiJ), 


;N    UA7A 


Ht/SC  (2f  ii9  J      U,  ,Y 

C'J     121     i  =  1  ,  i\  X 

R  t  MJ  (  ^  ,  1  -  J  )     \  A  (  1  ) 

121 

CCM  IN  II 

0")       12*      I=l,vlY 

(<E:i;i2,ijjj    ya(  i  j 

122 

CCNTlMc 

00     l^j     I  =  1  iNX 

CO     12  3     J  =  i  ,   !  r 

S  E  £  0  1 2  ,  I  2  0  )     ^AIJ   ,J 

123 

CO  J!  INu-8: 

CO    221     1=  L,   IX 

R  E  a  0  (  ?  ,  I  3  u  )     X  o  U  J 

22  1 

COOT  IN  J  6 

l)C    2  2c     i  =  1  t  N  Y 

RE  AO  (  2  ,  1301     Y3(  1  J 

222 

CC\T^  JE 

OU    22  3     i  =  l  ,  NX 

Ofi    223    J  =1,0  Y 

PE*C  12,  12j  )    I'M  1,  J  J 

223 

CC  ^  TIULfc 

C    T  R  Al 

£0    2  2**     i  =  1  ,Nf 

224 

W>II|sUlil 

CT    <^b     1=  1,   '* 

221 

YAM  i)=X/!  1  1  ) 

0  1     2  2c    1  =  1,   J  Y 

00    IJ.I     J=i,  !< 

22t 

ZA>(  I  ,w)=Za(  w»i  J 

00     Zc7     1=1,   !Y 

227 

XbXi 1 )  =  YE  l  i  i 

00     2  *b    I  =  1  ,NX 

228 

Y^X{  I  J=Xfc  (  1  J 

CU     2  29      1-1  ,  NY 

00     2<S    J=  1,  -1  a 

229 

2uX  (  1,  J  1  =  /  .'.  (.1,  l  ) 

C    F* 

CBt     ShlliNui! 

YP  A  =  C  .0 

YPl'=bb. 

XPA=C. 0 

XPC=2C  . 

XP0  =  -2  J  . 

C     PR 

OGSAf    CONSTANTS 

YrLC=-jO. 

YRIF=30 . 

YRI  M=-3C. 

XPLC=-30  . 

>RIP=3  J. 

XR  IN  =-30. 

ICCFS=1 

.NiJCuFS  =  2 

NOSl*=  1 

ii  r  =  i 

ISCAN=1 

PEL XP5=  J.  3 

C    tXP 

19j: 

C  0  0  v  r  '  L  E 

401 

FLPN^T  (r15  _/t  j 

10  7  J  JO  10 
INTG0020 
IN  TClC  iJ 
1  IT  J  jO',  J 
I N  T  C  C  C  5  C 
INTOjO'jO 
I  0  T  0  J  J  7  J 
INTcOOfiG 
I N  T  C  C  0  <j  J 
INtGj1>jJ 
INTlCUO 
IN  TOO  120 
I OT0J1 JO 
10  1  CCl'iG 
INT  0  01  30 
10  TCClcC 
INT0J1  /O 
I  N  T  J  j  1  ) 
IN  Tu  ISO 

int  00200 

INT0U210 
I  N  f  0  0  *  2  J 
INTO 02  3  0 
INT3J240 
INTCL2  5G 
I  ;<T'jj  2dJ 
I0TCC2  1C 
INT  002  dJ 
IN  TO 02 90 
INTCC3C0 
INTOGilJ* 
10  TC 0  3 20 
IN  TOO 3 30 
INT  00  3*0 
IN T CC ibC 
IMTCCitO 
IN TC 03  /O 
I N  T  L  C  3  e  0 
1\T  003  VJ 
INTCJ400 
INTCL41G 
INT00420 
1 N  T  C  L  4  3  0 
IN  TOG  140 
I  N  T  C  04  b  0 
IN  TC  CcG 
INTU04  70 
INT004oO 
INTU4SG 
IM  T  J0500 
IN  TO:>b  10 
INTCG520 
INTJJ5JQ 

INTCCbIC 
IN  TOO 5 50 
INTCJSaU 
INTC0570 
INT  00b  00 
IMT0ob90 
INT  liOcGG 
I NTOoo  lJ 
iNTCt£2C 
I N  T  0  C 1 3  0 
INT0J640 
I N  T  C  C  L  b  G 
IN  TOOooO 
INTJ06  70 
10  TCLcdO 
10  T Out. 9 J 
I  N  T  C  (.  7  C  0 
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FILE:    ^VJfi 


FORTRAN 


REAO(b,4JlJ     PT 
WR  IX t  (,*>,  -.u  1  )    -'T 
If  iNiI^  1  w  •  c'j  .2  ) 
REAl)(  b,  iCI  i     PS 
* R  I T  £  (  oi  '+ »  i  J    °i 
RtAU  (*>  , 4ui  i      ^Lh  A 
*f%  ITfl  ~,  4ui  )     ALF/ 
FE/SC  <b  i-tJ  1  J     Pl-ii 
WR  I  7fc<  t>  »'->01)     FrU  i 
CALL     i  J  TPL  i  (AAf  f< 
AlFC=ALF£-XFC 
CALL     I   J  IVL  I  <  aA,  r  : 
ALFC=AL  M-A^„ 

Call   iintflT  (XA»y.' 

P  M  I  P  =  P  H  I  i  -  i  ?  : 
CALL     IMTCL"!"  U"iY  t 
*RI  7L  (6  ,'+-J<l)     CP-  , 


L    TC    5 JO 


,H  t  ALFA  ,PHI  I  iNXtNYtCBA) 
tZA  ,ALFC  f PHII  ,\X,NY  ,CPC) 
,ZA  ,ALFD,Phl  I  ,.\X  ,  NY  ,  CRD) 


402 

(•.KMfTI  •     T-  Ub    LPi 
FJYIn  =  PT-FS 
PPA=CPA*FjY\+P5 

fp^cp  e*°c»  '  *-°  5 

FFC=CF„*  Fjf  f.tPS 
t>PD=CP  J*P^Y"|f'"i 

I'l^ir  lu.  -j\ 

WRi  Tt  (o  ,4  )  1  i     FrA  , 

FFE,FPC, PPD 

402 

FJRMA1  (  •      il  'ILL    r*F 

S  •  ,  4  F  1 0  .  ?  ) 

C    RE£C    ME^SURt^E     I'i    OAi 

'  A 

IF  (NCSiV.L',.1  J     o. 

:  to  sol 

50C 

READliil3SJ     ">da 
WRITE(o,i_.o)    r,-"\ 
Rt;C  lb  ,135  J     ?^ 
WR  I  rt<  u,  1  Jul     PPH 
R  1 1  £  I  b  ,  i.  J ')  ;     F  F  C 
WRI  IE  lo  ,  1  J  (.J      Pr^C 
REiClb,  1  ij>  )     'PJ 
WRITE  (o  f  lio)     PB0 

£01 

REAJ  I  tt  1  3  5  >     P  ii  -J 

wRIT  E  (o,  IjG  J    P  j  1' 

I F  (  I  5  C  A  J .  l  J .  1  J     o  - 

TO     181 

18C 

IF  {  I  SCAN  .G1.1C.J 

1IT=1 

ICCPS  =  i 

ISCAM=  ISCAm+1 

FS 1 N=  F3 I IW5.  J 

*Rl  TE  (  fc,  li tJ     PS  IN 

J>iJ    70    305 

13  1 

CCMIf  Jt 

FS^PSIN 

PSC=PS 

fst=ps 

C    CALC 

CEFFlCIfcNTS 

30  0 

co.<T  lr:  jf 

WRi  Tt t  j  t'i'iv  )      i  IT 

S4C 

FuRMATl  ' C  , '     1  II - 
FUY  J=PT-p J 
CPA=(FFd-hi)/PJYi\ 

CP3=(Pp  i-'j]/Pj  r  . 
CFC=  (F^'.'-PS  )/^>  „i    : 
CPO-l  -PJ-FS)  /rJY.". 

ATI C,\    0j.     • 

900 

^irnJf  wj    c^ , 

CM.f.PC  t'°0 

912 

FCRf*;!" (•     LLC  J    lPS 

'  f  4  f  1  J  .  f»  l 

C    CORR 

i\      T'_      E  '  i  jOr.C 

IF  (NCC  JPS.lO.^J     l 

i  "n   jji 

2AI*AX  =  -1C000. 

.CALL    ''l/r.UZJ.Z^ 

'\  a  ,  JX  t\'Y  ) 

IFKOJ.CT.Jl.AKj 

3  7C-- 1  .0  i.*P3 

IF  (CPC  i  I.  i  i  >!:,<) 

f  ; ;  - 1 . )  i  ■■  f  s 

1F{  CPU.  J  i'.,:  i   'A  \  J 

PjC^I.JI'-PS 

ZUVAX  =  -1-J^j  J  . 

CALL     Vu^MZ'i.Zc 

h-iX  |NXfi!Y  ) 

IF(C°fc.CT./i  iaa  ) 

1  SC=  i.  .0i*Pi 

ZAH  r=l03J0. 

i.Zf.i  "-LFA  ,  3HH  ,NA,NY,CP*) 

C  F  U  i  C  F  L  ,  L  F  0 


.15) 


Lr 


b     IN    CAl  MR  AT  I 'IN     RANjE 


Tj  fa  J  /  U 

1NTCJ7/U 

Iiyl  TCL7  jj 

1NTOL74J 

I  NT  Co  7  "30 

IN TC07tO 

l.'JT  Ju7  7J 

I'J  TCl  ?e  J 

I'll  ll/So 

rujj?oj 

I M  1  L  C  6  1 0 

IN  V00 e<  20 

I  N  T  0  L  fc<  j  J 

IN  7  Cue  4 J 

I N TOO b 30 

IN  i  C  u  o  uO 

IN  TOui 7 J 

I.NTOJot'.o 

i  N  7  C  L  d  S  J 

INTuu^  JO 

IN7C  JOiO 

I N  T  C  C  9  2  C 

ImT  00-/30 

H  T  L  o  S  4  J 

IMT0u9  3J 

I. NT  009  60 

I  N  7  C  C  c.  1 J 

INT  009  dO 

IN7CGf;90 

IN  TClu CO 

INTO  10  lu 

I  jTLIC'u 

IN  TO  10  10 

INTClo  '♦  J 

IN  TO  IC 50 

IN  TO  10  60 

I  rj  to  10  /o 

INTO  10 40 

INTOlu'.O 

IN  TO  11  CO 

INTClilC 

INTO!  1  A) 

INTCIHU 

1NT011A0 

INT01lt»L 

IN  7C  1  IrC 

INTOil/j 

iNTOlidJ 

IN  TO  11  >y 

1  NT  )l<_00 

IN 7C121C 

INTol220 

intu<;  30 

IN  7C1^hU 

INT  )L^5J 

IN  70  12a  J 

I  J  ITC12  7C 

I  NTO  1*.  oO 

I J7C12SC 

INTOi-jJ 

INT  oi  J  10 

INTL1J20 

INTO  13  30 

I  NT  J  I  jw 

I N  f  C  1 1  5  0 

I  NT  0 1  3  6  J 

I N  T  L  1  J  ?  0 

IN  T  CI  3c  0 

I NTui j9j 

INTCl^CO 
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FILE:    v»VCF 


FORTRAN       PI 


'i»*i' 


963 

96  4 


S2L 

921 

301 
C 


,   15) 

TC    180 


•,F^0.5) 


CALL    MHA^YIZA.ZAMN^X.NYI 
IFICPft.LT.Zaw4)    P$C  =  o. 99*(>S 
IF  (CPC.LT.  <.AW  M)     FSC  =  0.99*PS 
IFtCFO.L  I.  ZA  Uw)     P$C=C.99*PS 

ze.WN*  u  jjj  . 

CALL    MNJfY  (Ze  ,ZE*iM  iNX  ,  NY  > 

I  F(  C-^  .L  r  .ZOMIN)     P-3C  =  C.99*PS 

IF  (PSC.EU.P".  T  )     uC    TC    3U1 

ICCPi=lCCPS«-l 

rfRITEIo,  »fc«)     ICTPS 

FCnHsTC      ICCPS     aS« 

IF  UCCPi.GT.  1C)     oC 

FS=PSC 

FS  7=  f-SC 

WRIU(o,S2lJ    PSC 

FOP'^TC    ^SC     IS 

GU     10    3  00 

CCM  INUE 
CALCULATES    T  W  3    ALF    fl     Pnl     ANGLES     FH°     •!'     6     Mi 

CALL     1M  I  cC$(  YPa  ,  YP-1  i*RlJ  lYRUPf  YXI  Ni 
lCP*»CPdf/»Lf»PhiiAafY/ltZA| 
2>B  .VHiZfif XAX,  YAX  w /X 

ALF  1=alF(  1  J 

£LF2=£LF  (2  J 

PHI  1=P,<I  1  I) 

°HU=Pril  (  <:) 
CALCULATES    TWC    -LF    I     Fni 

CALL     1  nI  T  5C  3  (  A P  ft  ,;<  PC  ,>  R  LQ  »  X  f<  U P  »  X  *I  N  » 
lCF^,CPC,Pf-i,  .iLF,AiX,Y4X,ZiXf 
2>*>»>AXfZ^X,XA,YM,ZA,X4fYA,ZA|NY|MX 

ALF3=ALF  1 1J 

/Lf4  =  ALf  (<d  i 

PH13  =  PHI  (  1J 

FF14=PHI  (2) 
SELECTS    The     FROPtR    ALF    L    Phi     UIGLE3    HUT 

FALF  =  (  (ALFi>  H_F2^LF3+4LF4)-AMAX1  {  A.LF1  ,ALF^!  ,ALF3, 
lfLF4J-JMIl(iLFi,ALFi,ALFjf.ALF't))/^ 

FPHI  =  (4P.-iil*-Fni2*Knij*P:-iIt)-AMAXl(Prili,Phl2tPHl3tPH^) 
1-AMINK  PHI  l.Pril «:,?«!  I  2,  PHI  4)  I  /2 

'^RlTEIo  ,  ;Oi)FALF,FJhi 

FJRNATI'     SLLtCTtu'  ,2F10.2) 
C    CALCULATES    .'lErf    Pi    F«UH    DATA    OF    POSITION    D    AND    CORRECTS    PT 

AlFlJt^rdLF-XPJ 

CALL     IiMTFLT(  AA,  f,>,ZA  ,A  LFO  «J  ,F  PH  1 ,  NX  ,  NY  tC  PD) 

ALFCN=F«LF-X°C 

CALL     INTFlT(XA,Y/!,ZA,ALFC  J  ,F  PF 1 ,  NX  ,  NY , C PC) 

WRITFio.SlO)     CPCtCPO 

FCF*<£T('     N  EV    PC     6    C^C    CALLEO ',cF15.  S) 

PSNMPl5C*CFO-PPO*CPC)/  (LFJ-CFC  J 

HR  ITF(t>  ,  *1  i  )     rAuJ   fF  PHI  ,PS'i  ,P  T 

FCFC;t('     ITl-.ATiCN    VALJcS'  ,2F6  .2,  2FL0.2) 
C    CONVEX  OEM  t      iE  i!3 
950         CtLPS=PS^-PS 
S51  WKl  Tt(  j,952  )     FS  ,  FS.'-I  iCElPS 

FORMAT!  •     EPSPS    V4  LUfc  5  »  ,3F  15.5  J 

EPSPS=  Ae  ">(  (i'Vl-^S  )/P  i  J 

l»R  I  TEl  o  ,-"j6l)     t  P3FS 

F0?M£Tl  •     tPS^S     I  .,  •     ,F25.15) 

fcPSPSb=.0UOUO000o099*(PT**2 
) 


,XcX  ,Y  EX  ,ZeX,NX  »KY  ) 


/HOLES     FOR    •  I  •     €    'III' 


; 


OF    THE    FOUR    VALUES 


900 
SCI 


SOS 
910 

90  7 
911 


952 

St>C 
961 


0 )-. 00000098 *PT 
IF  (ItPjl.fcO.^)     crJSPSG  =  .C  COO  GOO  COOS*  I FT#*2. 0 ) -. 0000  003*PT 
IF(ltP^.t').J)     LPiPiL:=.COOuObJLJj9*(  PT**Z.  0)  -.  CCoOv.CE*PT 
IFItFSPS.cT.EPSPSbJ    GC    TO    310 
-GO    TO    30  I 
310         CCI^TIMJE 
I  I  1  =  11 i  +  l 

IF  HIT. GT. 10)     v.0    T3     205 
FS=PS+RtLXF5*(P£N-PSj 

psr.=  0i 

PST=PS 
$24  *RiTt(6,925)     PS 


INTO  14  10 

I  N  T  J 1  <t  2  J 
IN TC  1430 
IN  l'01'*4o 
INT014  50 
IN  T  C  14ZC 
INT  0  14 70 
1 NTC148J 
I  MYCINS  J 
IN  TO  1500 
I N  T  C  1  5  1 G 
INT01P20 
I  NT  01  5  30 
IN  TC  1540 
INTO  13 50 
INTO  15  60 
INTU1W0 
I NT01560 
I  N  T  L  15  90 
1 N  T  0  1  o  0  J 
INT01610 
INTCU2C 
INT01030 
INTu  16  40 
INT  CUSO 
INT01660 
I N  T  C  1  t  1 C 
INTOloSO 
INTC1690 
1MTCI  no 
IMT017  10 
INT  CI  720 
INTG1720 
INT01740 
INTC17*u 
INT01760 
IN  TOW  70 
INTC1780 
INTO  17  90 
INTC1EC0 
lNTGiaiO 
INTO  Id  ^0 
INTC1830 
INTO  18  40 
INT01850 
IN  TC IE tO 
INTUL870 
IMTClbdO 
INT0H90 
I  NT  01 900 
IN TC  IS  10 
INT  CI  920 
INTO  19  30 
INTC1S4C 
INTO  19 50 
INTO  1930 
INTO  19  70 
INTO  19  30 
INTC1SS0 
INTO2OG0 
INTu^OlO 
INTC2C20 
INT02030 
IN  TC2C-*0 
IN  TO  20  50 
1NT0^060 
IN  TC2C  /C 
INT  020  30 
INT02090 
INT02100 


74 


FlLtl     V>VGh 


FCRTHAN       Pi 


N 


S2  5 

303 
171 

3C5 

30/ 


15C 

C    CH 
C    CA 


F3RMATM     CURn.    PS    Ii    •  ,Fl  5 .  5 )   .u»i- 

GO   *C    iJJ 

hrlTtb,WD     F«aF,FPH,P3N,PT 

F  JR-1A  l(  JF13.2.F  15.3) 

GJ    TO     19  3 

WR(  TE(t>  ,  10  7J 

FORMAT!'    iC^HEO     100    TIMES') 

G3    TC    Hi 

E  NO 

SU*R0lJTlNc     I  lTJCi(iP.Pl,4RP2|4.*LTf  ARUPtA^INi 
1HT1  tHT2tt«eSi»PtS*»>U»Ylf  21  . 
2>2,t^.Z2fX^fYJ,iiJ.AttY4,Z^.M,N^) 

UlVfcNSiCN    :l5  1(  lu  ),"-  Ei2(  IOJ.'.RkE  >(  10)  ,  a!(4J)  ,Y1(  <td)  , 
l^l  iU',4C)  ,  X2  (<»  Jl  ,Y2  140  J  f  22  14  J  ,v.)  )  ,  X  3  (  4  J  )  ,  Y  j(  40  ) , 
2Z3(4C,<C).ai(-»0)  , YV< iC)  .Zt (40  ,40) 

PtSi  (1  )  =  -iOO  J. 

RE  <1 (2)  =90o0. 

R  E  3  2  ( 1  j  =  -J  j  J  o  . 

RcS2  (2)  =9JJJ. 

ISL=1 

1SR=  1 

0AR=1. 

AR2R=ARI  N 

GG    TC    Vj<. 

4R2aj=AR2< 

AAM  =  AP2F-A«iPl 

ARv,2  =  Anit^-;Kr'2 
ECK    FC^    LCrfeR    C  «L  I 

IF  UriNq  ,LT.^L0) 

IFIA3M2.LT.  V-L.J) 

K=l 
LClLATtS     INITIAL    PuNcTf'ATlCN 

CAH.    PfcUP  I'M   UfN. 

PAJL  =  £.<Rt>  (1) 

RA JK=ARRE S(2) 

CALL    Pt'JPT.i  (  J  1,1  2,X2,Y2,Z2,!-IT2,AR-12,ARRES) 


FRATI3  4    RANGE 
GJ    TC    150 
GJ    TJ    15J 

PCI  f*T.S 
,Al,Yl,*.l,rlTl,4RMl,ARRES) 


RB«L  =  Ar<?l:3  (1  ) 

RBJR-ARP.ESt  2) 
C    CALCULATtS    SUCctSlVE    P  EN  ETRAT  1 ON    £>CINTS 

DO    140    i-1  ,10000 

AR*R=A  *2R04R 

*R*i=A«2  rJ-ARPi 

ARM2  =  /i"<2R-A  <?2 
C    CHECK     F3R    uPi»Ei<    CALIBRATION    RAN^E 

IF(*P"l.^T.  J-UP)     GJ     TC    195 

IF|ARr'2.-JT.MKUP)     ul    TC    19a 

UiXl,  VI,  2  1,HTI,AR>11  ,ARRES) 


C    CH 
164 


CALL    Pt'JPTSru, 
RAI L  =  AShES  li ) 
RAlR  =  A.<ie.S(  2) 
CALL    Pt^lPfS  HI,  ')2,X2,Y2,Z 
PdiL=*&ScS(l  I 
RBIR=ARRES( 2) 
IF  (  ^AJL  .E  J  .i'JJJ  .  )     ISl=  i.0 
IF  (Rt JL.t  ..1000.  J     l3L  =  10 
IKRPJR  .L-J  .  IJJJ.  )     13R=10 
IF  I  R£  JR  .  c>.  .  1 000  .  i     I3P=lu 
IFIRmIl.E...  1006.  )     1 3L  =  1J 
IF(RPIL .EO.luOJ.  i     I SL=  10 
IF  IRPlh.  tC.10  >0.  )     1SR=1U 
IF  (?AI' .  E'J.  1000.  )     ij':  =  lO 
fcCK     FCP    LlUCT     FIT     J  I     LEFT 
IFi  1  SL.  NL.10)     GC     T3     lo4 
-ISL=  1 
CL    TC    18J 

IF  (  «M  l.  (Mt=  -  RBI  L)     GO    TC    180 
RESKis  )=R'UL 
t>kc,2  (K)  =  *k«iR 
IFI*.E«.2)    Gv    TJ     155 
K=2 


2,HT2,iR42,ARRES) 


INTERSECTION,  EVALUATE    IF    RE011RE0 


IN  TOil  10 

INTQ^I^O 

1N7C2  1K 
1  JV02l<tO 
I/MTJtlbO 
INTC2  UO 
INT 02 170 
IN  T  J  2  1  £  G 
INT C2  ISO 
INT02200 
1 N  T  C  2  2  1  0 
INTO  22*0 
INTU22 JJ 
IN7C2<:40 
1MT022'JJ 
INT022oO 
INT022  70 
I  "-IT  022  JO 
I N  T  C  2  2  S  0 
INTC2iC0 
HT02J10 
INTC2J20 
INTC2:>30 
INT0234 J 
INTC2  3SO 
INT02J63 
INTC23  7C 
IN  l"0  2  JbO 
INTC2J90 
INTC24C0 
1NTJ2410 
INTC2420 
INTC2<«iO 
INT j2h40 
I N  T  C  2  4  5  0 
INTC2*»eO 
INT024 to 
IN TC24EC 
TITJ2h  }J 
INTJ2  5 J T 
INT  02  5  10 
INT025  20 
INT0253J 
INTQ2  54J 
INT0255J 
INT'J25->J 

rirc257o 

IN"f02SriO 

INTC25S0 
IN  TO  26 00 
INT02(«1J 
IN  TC2(20 
INT02O30 
IN  102640 
INT0265O 
INT  02  6 oO 
INTC26/J 
IMTC26E0 
INT02O90 
INTG2  7CO 
INT  Od  HO 
INTC2  720 
INTC2  110 
I  TJ  T  itlfnO 
INTC2  7  50 
INT.J2/60 
INT  02/  70 
INTC27aO 
INTC27S0 
IMT02B00 
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FILE:    WVD° 


FORTRAN       PI 


GO     70 
ChECK    (Tk 


Id; 

i  UERMEJI  4TE 


•»!*•  ■ 


GO    TJ    160 


180  IF4lSL.Nt.lOj 

ISL=  1 

CO    TO    idi 
16  C  IF  M  F>-JL--.eJL)  /(  I-  A  iL-Rrt  I  L)  .GT 

AA=(  AF  2RJ~  W2R  )/(  R;jL-r.AlLJ 

EA  =  AP*R-*/>*RAlL 

A3-(  AK2SJ-At<2Kj  /(FdJL-KSlLJ 

B3=AR2R-AB*FFUL 

KESl  (KJ  =teb-FA)/  (*A-*E  ) 

REG2<k)  =  AA*  <t  jU  is)  fflA 

IF  (K.LC.%2)     G  I    TO     1;>5 

K  =  2 
C    CHECK    rJR    DIREC1    HIT 


LEFT    INTERSECTION, EVALUATE    IF    REQLIREJ 


0.01     GC    th     id3 


n 


XIGHT 
lol 


183  IF  (  I^F  .M.iO  J 

I  S*  =  l 

GO  TO  16  J 
161  IF  tK^lrf.lvfc.seiRJ     CC    TC    182 

KE  Hi  KJ  =  *H  R 

PES2  (K  J=A^9 

IFlN.fcJ.2j     uJ    TJ    155 

K=2 

GO  TO  16  J 
C    CrlECK    FOR    I  i\7E  RMfc  Jl  ATE 


IiMTEriSECTl  CN,EVALUATt    IF    RECUIRfcD 


18^ 


162 


IF  USF-.ME.  U)     G^     n     1 1 2 
IS  B  =  l 

GO    TO     lei 

IF  (  <R^J  R-Rb\J  <,)/  RiF-^nH  J  .GT 

4A=(*P2Ro--*cf<)/(FAj;<-KAIRj 

BA=',R2R-'U*->  AIR 

*B=  (/si-.^J-^^  J/  IR3JR-RB  l\  ) 

BtJ  =  AR2K-A3*R'iIR 

RESKK  J  =  ('l-*-4aj/(  A<\-aB) 

RLS2  (KJ  =  AAVrl:S  1(  K  J  *<z£ 


RloHT    INTERSECTiOMiEVALJATE     IF    REQUIRED 


CO)    GC    TO     163 


IF(K.E  y 
K=2 
RE-IM  TI  ATE 


l\ 


J    TJ     15  5 


lt>3  Wi3j=AR2A 

RAJL=FAIL 
R  A  wR  =  3  A  i  .-i 
RBJL=RB  IL 
Re^K  =  Fb  i? 

14  C         CON  UN  IE 

WR  ITEJu,  21  J) 

210         FCF^Tl'     LOOP 
GO     fO    155 

195  !F(RESi(lJ 

IF  iRcSL  (2) 
IF<RFS2(  I  J 
IFlPfcS^  (2) 
GO    TC    15  5 

C     EVALJATICr      )F 


140    TFSajurii  ND    PILOTS     FOUND1  J 


EO  .9  000.0  J 
c  J. 9  J 00.0  ) 
l  )  .  y  v>  w  J  •  J  ) 
t  J.9  0GJ.0  J 


JO 
J  3 

G  J 

a 'J 


TO 

T  w 

r  j 
r  ] 


191 


VGC=0.  J 

CAlL  PCNP'I 
Rl=;  tS  (  ifif 
CALL  H:.\FT 
R2=  A3S(  A3> 


INTERSECTIONS    WHE' 


19  1 

191 
191 

191 

CLRVES    ARE    ALMOST    TANGENT 
D,APRES J 


( N2 #N1  tX3  i YitZi  tHTl  »Y 

L  S  ( d.  J  -  h  ?  R  E  i  (  x.  J  J  /  2 

">  (i\2  iul  ,  <<+  jYh  ,  z.  t  ,n72  »Y  u  J,  AF  ^ES  J 

F  S(2J-~RRt  S(  1J  J  /2 
FES2  U  )=(  IK  l**<ij- (•<<.<**:)♦(  A:<P2**2J  ) /I  2.  0*'.*P 
RE  £2  (2  J  =RL^2(  1J 

RESH  1  )=SO\T(  AroU  Rl**2)-Ut  S2(  1)**2))  ) 
FES1 (2J  =-RLSi(l J 
15  S       -RETlxN 
ENL 

SU0RCU7  I  l\t     FlNPTS  I  MX  ,  NY  ,X 
Di^ENSI  JiM     yliil.LMiCliZI 
IR  =  1 
C     SEARCH    FCF     o    LF     u^Ei-    V     LINE 
DO     lul     1=1, OY 
IF     4  Y  (  1  )  .U.Yu)     oC    TC     lo^ 


,Y  ,Z,  FT,YG,XR  ) 
4i,  ,  40)  ,X(  40  J  ,  XI 


(1  JJ 


NT  J*  J  10 

NTo.iH20 

n  t  o  i  a  j  a 

NTj^j'.O 
N  T  0 2  b  dO 
NT028fcJ 
NTo2b  70 
NTC2ccO 
N  T  C  2  E  S  0 
NT 020 00 
N  r  C2 S 1 0 
NT02920 
NT  029  J  J 
NTC2S40 
NT0295U 
NT  C 2 c.  1 0 
N  T  0  2^  7u 

nto^oso 

NTC2SSG 

'■J  T  0  3  0  C  0 
NT03J10 
NTG3G<:0 
NTui030 
NT0j04J 
N TCJC5J 
NT0J060 
N  T  C  3  C  i  C 
NT  J3u80 
NIC J  COO 
NTC31C0 
NT031 1J 
NTG312J 
NT0313Q 
NT0J14J 
NTC3150 
N  T  0  3  1  6  J 
NT.jJl  70 
N  T  C  3  1 1 C 
NTUJ190 
NTC^2CJ 
NT C 32  10 
NTQJ220 
U  T  C  J  2  3  0 
N  TO  3  2  40 
NTL.J250 
N  T  C  J  2 1 0 
NTu.i27J 
0  TO  32  30 
NTuj290 
NT 0  33 00 
NTC33  10 
N  T  G  3  3  2  0 
N TO 3330 
.'  I  T  C  3  3  4  J 
N TO 3 3  50 
.N  TOJioJ 
H TC j3  70 
NT0333J 
N  T  0  J  3  'O 
NT0J400 
NT0J410 
'J  T  C  2  4  2  0 
N  T  C  J  4  3  0 
NT0j440 
NTC3450 
N  I  0  j<tcJ 
N  TO  34  70 
M  T  C  3  4  6  0 
NT0j-»'9O 
NTC35CC 
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FILE:    WVufi               FORTRAN       PI  MP                                   S 

1C1         CuNTIhLE                                                       if»*iv  INY035  10 

102       j=j-i  imi^^o 

C    CQAkSL    SCAN    FCh    ZERO    P*SS  •                                               1NTGJ520 

C    NEXT    <♦    LINES    Vt    EXfcClfEG    IN    FIRST    SEARCH    CNLY  1NT)J5<0 

YC=(YC-Y(Jj)/<r<J-»x)-t{j))  I N 7^3 3 j J 

ZXil)  =  UlliJ4lJ-..Ui..)J*YOZ«ltJJ  l.lTUJSuO 

IFCZ/(  i  )  .fc'J.ifl)    ,i    7:     U/  INTOJ57J 

Id  =  ^  INTCJSfcO 

112  JJ  1  C3  1=1-,  ,  i\X  rJTC35S0 
ZX(  i  J=(  Z(  1  ,J  +  1)-Z(  l»  J)  )*YC,4Zl  I  tJ)  I.NTOioOO 
ZM=U(l».,«-x)-Zll-i.,Jji*<YJ-Y(J)J/(Y{J*l)-Y(J))+ZU-l.J)  INTC3MC 
IF  IZ<(  1  )  .  c'J.H  I")  oC  l'C  108  INT0Jfc2G 
IF  I  (/!Xl  I-U-hT  i/  (ZX(  1)-FT)  .LT.c.G)    GC    TO     10<t  iN7u3c>30 

103         CJNTlN)c  INTCJfc^C 

GO    TO    lit  lNT0i6b0 

C    CALCULATION     TF    DEMETRVUDNS     I  J    EITHER    TRIANGLE  IN7G3600 

1C4          *M=M1-1  )*lxi  x  )-a  (i-1  )  J*(Yo-Y(  Ji  J/  (Y  U*i  J-YU  J  J                                                         IMTCJt/J 

1F(Z'4  .dO  .h  f  )    G'J     "i  3     Kj  IMJJooJ 

IFU/xm-H7j/Us-H|.LT.J.U)G3TO  10  3                                                                               INTC36S0 

DO=A^S(  n-^.V)+Abr  (ZX  (x-1  )-nTj  INT  L  3  7  CO 

XS-=X(I-1  J*Ia^-a1  I-IJ  )*io  j(Za{  I-l)-rtT)/O0  IN 70371 J 

GC    TC    110  HTC372C 

1C5          JS^X1+(  x(I  )-XM)*A3,»(  Z*    -HT)/(At:SU"l-HT)+AB:>UX  U)-nT)  J  1NT037jO 

GO    TC    llu  IMTo.i7'tO 

C     LOGICAL     Fti\fPATlu.NS     ACCL  f  ULAT  1  CN  INTCJ75U 

10  i         XR  (  1  )=  >>(  1 J  INT.)37tJ 

1  =  1  1  N  T  C  3 / 7  o 

GO     70    113  INT037HO 

10  8         XR  (  P  )=X(  I  )  I  NT  03  7 -iO 

bu    tc    111  INTC38C0 

I'll          XR(I3)=><  INTUjoIO 

111          IFUP.tO..)    GO    TO    114  INTu38^0 

113  IR=2  IN  703830 
lR^I-tl  1NT03c40 
GO    TC    112  :       lNTOiBbJ 

lit          IF(IR.EC.I)    oC    TC    117  lN'lC3ttC 

XR(21=10oj.U  lUTOitJJ 

GO    TC    ill  IN7C3tfcC 

111          >R(1)  =  1COO.O  INT03890 

XR(^J=1C0J.J  '                                                                     INT03900 

114  PETUPN  INTC3S10 
ENO  I.NTC3^20 
SUeROUTlNE     IITPlTIX,  Y,Z,  X3,  yG,Ka,  !YiZR=S)  INTGj^iO 

C     LINE  AW     ^TtHFCLATjCN     CI     CALIF*  AT  ION    SURFACE     TO     EVALUATE  INTC3S4C 

C    CP    VpLLE    A!'    CJ-QHOlN«TcS    >G,Yj.     Tn  b    RESULT    SETUKNtJ    IS    ZrES  INlOiSbO 

OIvENSx.rj    X  (tOI  ,Y  (-*J  JtZItJf  40)  INT03VtO 

C    SEARCH    FCR    I     LF     LEFT    X    LINE  INTC3S7G 

DO     1     IC=l,'iX  INV0J98J 

IFIXilO.CL.XG)    uC    TC    Z  Ii\iTC3St-C 

1  CONTINLb  INTo^OJO 

2  J=  iC  -1  INTOhOLO 
C    SEARCH    FCR    J    Cr     LCWEF    Y    LINE  INTC4C20 

DO     3    j;=l,!Y  IN T040 30 

IF  (Y  UC)  .cC.fGJ     CO    TO    4  INTG4U40 

3  CCMlMc  INTG^JbO 
*                J=JC-1  INT  04 060 

AL=  (YU-H)-Y  (J  ))/  U(  Hl)-X(  I  )  )  INTo4G70 

BL= V< J)-AL*XlIJ  IMTC<tC8G 

tCR= (AC *XG  )  +3L  '                                  INTU40yj 

IFjYG.cUa?)    u]    Tl     /  INTC-UCG 

IF  (  YG.i  1.  YC  U     GJ    T3    5  INTO*  110 

C    RESULT  It.Q    Z     IS     I  I    LJ*RER    TRIAJGLE  INTC-+L20 

X3=Xd  +  l)  IM7C4130 

-*3=t(J)  INTO'titu 

2j>=Z  I  i  +  l  ,  J  )  lNTj^l  5J 

GJ     TO    6  INTO-'t  loJ 

C    RASULTlIu    Z      IS     IN    UPFEK    TRIA.Nl.LE  1NTU-U70 

5                X3=X(11  IN  TO*  180 

*3=Y(J+1)  INTl^ISO 

Zj=Z11,J^1)  /                             INTJV20J 
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FILE:    WVC>P 


FORTRAN       PI 


C    GENERAL    CAlC ILATI JN SGC JU    FCR    oCTH    TPIANGLES 

6  X1  =  X(IJ 
Y1=Y ( JN 
Zl=Z(  I,  J  I 
X2  =  X (  1*1  i 
*2  =  Y<J*U 
l<>=Zi  1*1.  i*U 

AP=U*i-<JiMYi-Y^)-(Zl-Z2  J  *(  Y1-Y3  »)/UXl-XJ)*(Yl-Y2)-(Xi-X2) 
l*i  Vl-YJJ  ) 

ep=  ( (zi-z«  i-*p*(  x  i-xi )  i/(  Yi-y<:i 

CP=Zl-AF*<l-~P*Y  i 
ZRES=AP*  /(J+8P*  iij*Ca 
GO    Tn    1JJ 
C    RESULTING    L    IS    CN    JIVUIHG    LINE 

7  Z*F3=(Zl  I*  l(j*l)-Z(  1  ,  J  )  i  *i  Xv»-X(I  ))/(X(Itl)-X(i)J*-Z(I,j) 

ioo       feting 

Ei'U 

SO* ° CUT  I  ME     MlxA^YlZiZ'HAiJX.JYj 
C    EVALIAT1CN    OF     N^Xl^U^     J  F     CF     CALIBRATION    SURFACE     IbT'U 

OIVtNilJN     Z(4G,4U 

CO    1     I=1|NX 

DO     1     J^l  ,NY 
I  IF(Z(  I,  J  i  .Jt.ZNAX)     Z'UX  =  Zl  1  ,  J) 

tfETURK 

ENJ 

SUB3CUT  I  .IE    MIN4°.Y(Z,Z'-1IN,NX,NY) 
C    EVALUATION    Cr     ("IMVJ*    CF    CF    CAL  I  *T<  AT  1DN 

DIMENSION    ZIh,HI 

00    1     1=1, NX 

CO     1     J  =  1  ,  f\  Y 
1  IF(Z(  I,  J  I.LE.ZMlN)     Zf<IN=Z(I,J) 

RETURN 

END 
C    END   OF    WVD* 


SJ"F4CE    MATRIX 


INTC42  1J 

I  NTu't2  2J 

INTCW  20 

INT 04* 40 

1NT0>2vjJ 

I N  T  U '« ^  u  U 

ItfTO'tl  /J 

riTCt2cJ 

INYC42S0 

INTCooJ 

I N  T  L43  10 

MT')<o20 

IN T 04330 

IwTG424.J 

INT0-O5J 

I N  "l  L  4  3  fr  w 

IN  7  j si  /j 

I  NT  O'ficij 

IN  TG^jSO 

INT J V4UJ 

I N  T  C  '-.  4  1  0 

INT  J4'.2u 

INT0443J 

INTG44  4U 

1NTC44  5J 

INTC't<i6j 

IN TCs4  70 

INTJt4dJ 

INTC44C*'3 

INTC^jCO 

INTO 4 5  10 

INTUs520 

INT0-ib3J 

I  NT  04  5> 4  J 

INT  C4  5  50 
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